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ABSTRACT 

A  software  package  was  developed  to  perform  conceptual 
design  of  submarines,  using  the  Computervision  CGP-20OX 
Designer  System,  a  turn-key  computer  aided  design 
hardware  and  graphics  software  system.  The  philosophy 
behind  the  software  package  is  to  keep  all  major  design 
decisions  under  the  control  of  the  design  engineer, 
rather  than  embedding  key  decisions  in  the  program 
a Igor  i  t hms. 

Modules  a.re  provided  for  calculating  weight  estimates, 
principal  characteristics  and  envelope  geometry, 
resistance,  weight  and  moment  balance,  and  the 
equilibrium  polygon.  The  package  interfaces  with  a 
pressure  hull  design  module  developed  separately  in  an 
O.E.  thesis  by  Marvin  Meade.  Interactive  graphics  are 
used  where  appropriate. 

The  software  requires  a  knowledgeable  naval  architect  as 
the  user,  but  does  not  require  extensive  knowledge  of 
computers  or  computer  aided  design  systems. 
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INTRODUCT I ON 


GENERAL  DISCUSSION 

Computer  Aided  Design  (CAD)  has  become  an  important 
tool  in  many  industries  over  the  past  ten  years.  The 
capabilities  of  CAD  systems  are  growing  at  a  rate  which 
defies  efforts  to  catalog  them.  In  the  aerospace  and 
automotive  industries,  as  well  as  in  large  architectural 
firms,  CAD  has  become  the  standard  mode  of  design,. 

In  the  area  of  naval  ship  design,  CAD  has  been 
adopted  more  slowly,  for  a  number  of  reasons.  The 
comparatively  low  production  rates  bring  the  issue  of 
cost  effectiveness  into  greater  question  for  ship  design 
agencies.  The  tremendous  complexity  and  scale  of  naval 
combatants  makes  the  design  of  a  flexible  and 
comprehensive  package  of  software  a  formidable  task. 

Despite  these  problems,  a  number  of  packages  for 
conceptual  and  preliminary  design  of  surface  combatants 
have  been  implemented  in  both  government  and  private 
sector   design   organizations.   The   packages   used    by 
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government  agencies  have  tended  to  be  batch-oriented 
ship  synthesis  programs  with  little  or  no  graphic 
output,  and  virtually  no  real-time  interactive  graphic 
design  capability  (the  notable  exception  to  this  is  the 
use  of  interactive  interior  ship  layout  software). 

This  lack  of  interactive  graphics  has  been 
partially  hardware  driven;  the  concern  for 
standardization  and  the  lengthy  process  of  government 
certification  and  procurement  cause  a  substantial  delay 
between  the  initial  availability  and  subsequent 
installation  of  new  hardware.  In  addition,  the  CAD 
process  is  foreign  to  many  of  the  prominent  designers  in 
the  very  tradition-oriented  world  of  naval  architecture. 
Many  design  procedures  in  naval  architecture  are 
somewhat  subjective  and  difficult  to  quantify. 

A  further  impediment  to  the  use  of  CAD  has  been  the 
tendency  of  many  early  software  packages  to  usurp 
designer  perogative  by  having  key  decisions  in  the 
design  process  "hard-wired"  into  the  program  algorithms. 
Such  software  does  not  utilize  the  expertise  of  the 
designer,  and  frequently  precludes  trading  off  design 
parameters  to  achieve  optimal  designs. 

In  the  area  of  submarine  design,  the  above  problems 
are  further  accentuated.  The  volume  of  new  submarine 
designs  is  even  lower  than  that  of  surface  combatants. 
Accuracy  becomes  critical;  where  feet  and  tons  may  be 
acceptable  units  for  tolerance   in   large   surface   ship 
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designs,  inches  and  pounds  are  frequently  the  desired 
accuracy  level  in  submarine  design.  This  is  due  to  the 
hydrostatic  requirement  to  balance  both  surfaced  and 
submerged,  as  well  as  to  the  premium  attached  to 
internal  volume  and  deck  surface  area. 

Submarine  designers  are  fewer  in  number  than  their 
surface  counterparts,  and  they  often  have  many  design 
relationships  tied  to  subjective  criteria  (i.e., 
"designer's  eye")  developed  over  years  of  experience.  A 
program  which  does  not  allow  the  submarine  designer  to 
retain  this  subjective  design  flexibility  is  not  likely 
to  be  accepted  by  the  design  community. 

SUBMARINE  CAD:  A  HISTORICAL  PERSPECTIVE 

Although  Computer  Aided  Design  has  not  been  a 
primary  tool  in  submarine  design,  numerous  programs  have 
been  developed  to  accomplish  analytical  chores  that  are 
subsets  of  the  overall  design  process.  These  include 
propellor  design,  resistance  estimates,  hydrostatic 
calculations,  and  structural  design.  These  programs  have 
been  primarily  analytical  in  nature,  providing  no 
interactive  graphic  capability,  and  are  usually  written 
for  main-frame  computer  systems. 

Several  packages  have  been  written  to  achieve  a 
complete  iteration  of  the  design  process.  These  include 
the  CODESUB  program  CI].  It  provides   numerical,   rather 
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than  graphic,  output.  CODESUE  was  developed  by  the 
Center  for  Naval  Analyses  in  order  to  aid  in  projecting 
future  submarine  design  characteristics  for  both  our 
navy  and  its  potential  adversaries.  This  program  was  not 
intended  to  be  a  detailed  tool  for  the  naval  architect, 
and,  unfortunately,  has  many  key  design  parameters 
imbedded  in  the  source  code,  which  restricts  its  use  in 
trading  off  alternative  designs.  Other  packages  "graft" 
together  sections  of  existing  designs,  matching  the 
largest  hull  diameter,  to  provide  a  composite  design. 
ASSET  C2],  developed  by  Eoeing  Computer  Services,  while 
primarily  a  surface  ship  design  package,  may  be  expanded 
to  include  a  submarine  design  module  in  the  future. 
ASSET  provides  some  low-level  graphics  output,  but  is 
not  a  full,  real-time  interactive  graphics  package. 

Thus,  there  is  a  clear  opportunity  for  the 
development  of  a  submarine  design  package  incorporating 
the  use  of  the  full  capabilities  of  currently  available 
CAD  systems.  This  thesis  represents  an  exploration  of 
such  a  design  package. 

PHILOSOPHY  FOR  THE  THESIS  PROGRAM 

The  software  developed  in  this  thesis  accomplishes 
a  "first  pass"  through  the  submarine  design  process  at 
the  conceptual  design  level.  In  conjunction  with  the 
concurrent   Ocean   Engineer   thesis   written   by   Marvin 
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Meade,  the  package  includes  calculation  of  weight 
estimates,  envelope  and  pressure  hull  geometry, 
resistance  (speed  and  power),  basic  hydrostatics,  weight 
and  moment  balance,  and  determination  of  the  equilibrium 
polygon . 

The  primary  philosophy  of  this  package  is  to 
develop  it  for  a  user  who  is  a  competent  naval 
architect,  but  who  may  not  necessarily  be  an  expert  in 
the  use  of  computers.  Crucial  design  decisions  are  left 
to  the  user,  rather  than  being  embedded  in  the  program 
code.  Opportunities  are  provided  to  override  program 
algorithms  where  necessary  to  implement  the  desires  of 
the  individual  designer.  In  addition,  if  the  designer 
does  posess  some  knowledge  of  CAD  systems  or  computer 
programs,  it  is  easy  to  exit  the  package  at  appropriate 
points  and  tailor  the  design  to  provide  more  detail  or 
alternate  geometry  and  analytical  procedures. 
Interactive  graphics  are  used  to  provide  a  clear 
visualization  of  the  design  under  development.  All 
program  input  parameters  are  couched  in  standard  naval 
architecture  terms  to  facilitate  the  user's  interaction 
with  the  package.  The  title  chosen  for  the  software 
package  is  CADSUE  (Computer  Aided  Design  of  Submarines). 

Input  to  each  program  module  is  explained  in  the 
chapter  documenting  the  particular  module.  Samples  of 
program  output  are  provided  in  Appendix  I  to  this 
thesis. 
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HARDUJAREs 
THE  COMPUTERVISION  DESIGNER  SYSTEM 


GENERAL  DISCUSSION 

The  optimum  choice  of  hardware  for  implementing  the 
philosophy  of  this  thesis  was  one  of  the  several 
high-quality  "turnkey"  CAD  systems,  all  of  which  permit 
the  use  of  interactive  graphic  design  techniques. 
Graphics  packages  for  main-frame  systems  are  available 
from  several  prominent  suppliers,  but  are  a  compromise 
at  best,  since  the  processors  they  run  on  are  not 
optimized  for  graphics-intensive  applications. 

Suitable  systems  are  marketed  by  several  suppliers, 
including  Computervision ,  Applicon,  Bendix, 
McDonnell-Douglas  Automation  Division,  and  others.  For 
this  thesis,  Computervision  Corporation,  of  Eedford, 
Massachusetts,  offered  free  system  time,  instruction, 
and  other  technical  aid  to  develop  the  program. 
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THE  COMPUTERS  I  SI  ON  DESIGNER  SYSTFM 

The  system  which  Computervision  provided  time  and 
instruction  on  is  their  basic  "Designer  System",  built 
around  the  CGP-200X  processor.  The  heart  of  the  system 
consists  of  the  processor,  a  high-speed  tape  drive,  and 
one  or  two  300  megabyte  hard  disk  units.  The  CPU  is 
available  with  several  increments  of  random  access 
memory  (RAM).  The  particular  system  used  for  this  thesis 
was  configured  for  approximately  1.3  megabytes  of  core 
memory.  The  architecture  of  the  CGP-200X  is  16-bit. 

The  Designer  System  utilizes  two  distinct  operating 
systems.  Easic  file  management  and  system-level 
housekeeping  chores  are  handled  at  the  "OS"  level,  which 
is  similar  in  use  to  the  operating  systems  on  typical 
mini-computer  installations.  For  the  graphics  operating 
system,  Computervision  uses  the  " CADDS  4— X"  environment, 
which  is  tailored  for  graphics,  and  automatically 
interfaces  with  a  powerful  data  base  management  system. 
Graphics  commands  are  entered  in  a  simple  verb-noun 
syntax,  such  as  "INSERT  LINE"  or  "INSERT  SPLINE". 

Input/output  (I/O)  is  handled  by  several  devices. 
For  system  administration  and  text  file  input,  small 
alpha-numeric  terminals  are  provided.  For  graphics  I/O, 
high-resolution  color  or  monochrome  terminals  are 
provided  as  part  of  a  comprehensive   workstation,   which 
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mcludes  a  thermal  printer  for  working  quality  tent  awd 
graphics  printout.  Also  provided  at  each  workstation  ar& 
a  digitising  tablet  and  a  display  control  device.  This 
device  controls  background  and  graphics  intensities,  the 
number  of  lines  of  text  displayed  (4  or  24  lines),  and 
dynamic  control  of  the  display;  zooming,  rotation,  and 
scrolling.  Output  of  finished  quality  is  available  from 
line  printers,  color  pen  plotters,  and  black  and  white 
electrostatic  plotters.  Large  format  output  is  supported 
by  the  plotters. 

PROGRAMMING  LANGUAGES 

Several  languages  3.r&  provided  for  use  with  the 
Designer  System.  The  primary  language  is  Fortran-S,  a 
subset  of  ANSI  Fortran.  System  subroutines  are  provided 
for  interacting  with  the  graphics  and  data  base 
management  operating  systems.  Programs  entered  in 
Fortran-S  execute  with  the  greatest  speed,  and  allow 
creation  of  new  graphics  commands.  Because  of  limited 
CPU  space,  the  utilities  required  for  this  level  of 
programming,  which  ar&  often  provided  on  main-frame 
systems,  must  be  separately  programmed  for  each  module. 
Such  utilities  include  linking  and  loading  capabilities. 
Because  of  the  length  of  training  required  to  do 
competent  programming  at  this  level,  this  thesis  was 
programmed  in  an  alternate  language. 
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Two  macro  languages  are  provided  for  the  system. 
Both  allow  creation  of  execute  files  using  standard 
graphics  commands.  VARPR02  is  a  non-compilable  macro 
language  which  supports  Fortran-like  I/O  and 
computational  statements.  It  is  a  flexible  and  effective 
language,  but  executes  rather  slowly.  NEWVAR,  the  other 
macro  language  provided,  is  compilable,  and  therefore 
executes  far  more  rapidly.  The  I/O  for  NEWVAR  is 
cumbersome  and  difficult  to  format.  All  I/O  must  be 
handled  in  text  string  form,  and  converted  within  the 
program  to  numerical  values.  Despite  this  shortcoming, 
the  execution  speed  of  NEWVAR  led  to  its  adoption  for 
use  in  this  thesis.  Ultimately,  ''the  algorithms  in  this 
thesis  should  be  translated  to  Fortran-S  for  maximum 
speed  and  I/O  flexibility. 
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PRELIMINARY  UIEIGHT  ESTIMATES 


IMPORTANCE  OF  UIEIGHT  ESTIMATES 

The  algorithms  employed  in  this  thesis  assume  a 
weight-driven  submarine  design.  This  assumption  means 
that  the  envelope  size  calculated  for  the  submarine  mill 
contain  enough  volume  to  displace  the  same  amount  of 
water,  when  submerged,  as  the  weight  of  the  envelope 
displacement  for  the  submarine.  If  the  desired  design  is 
known  or  suspected  to  be  volume-limited,  an  arbitrary 
excess  of  lead  may  be  specified  during  the  weight 
calculations,  thus  driving  the  program  to  calculate  a 
larger  envelope  size.  The  designer  may  translate  this 
additional  weight  into  volume  manually  to  allow  for  a 
volume— 1 i mi  ted  design. 

Since  the  estimated  weight  data  will  drive  the 
basic  dimensions  of  the  submarine,  the  weight  estimating 
algorithms  must  be  made  as  accurate  as  possible.  The 
designer  must  carefully  consider  the  input  data  to  the 
weight  module,  and,  where   appropriate,   may   decide   to 
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override  the  programmed  algorithms  and  specify  the  value 

for  a  particular  weight  group.  The  program   is   designed 

to   automate   the   calculations,   and   not   the    design 

Judgement,  for  the  user.  A  knowledgeable  user  may  easily 

modify  the   algorithms   themselves   to   suit   particular 

needs. 

U1EIGHT  ESTIMATING  ALGORITHMS 

The  accounting  system  chosen  for  weight  estimation 
is  the  SWBS  (Ship  Work  Breakdown  Structure)  convention 
used  by  the  Naval  Sea  Systems  Command  for  ship  weight 
records  C3].  This  system  consists  of  seven  major 
numbered  groups  for  all  ships.  For  submarine  weight 
reports,  several  special  weight  categories  a.r&  added.  A 
brief  summary  of  the  weight  groups  is  as  follows: 

DESCRIPTION 

Structural 

Propulsion 

Electr  i  cal 

Combat  Systems 

Auxiliary  Systems 

Outfit  and  Furnishings 

Armament 

Sum  of  I -VI  I 

Total  lead 

Al  +  Lead 

(fuel,  stores,  etc, ) 


GROUP 

I 

II 

III 

IV 

V 

VI 

VII 

Al 

Lead 

A 

^ar  ia 

ble 

Loa 

d 
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NSC  Normal  Surface  Condition 

Reserve  Buoyancy    Main  Ballast 

Submerged  Disp.     Self  explanatory 

Free  Flood  Non-buoyant  flooded  volume 

Envelope  Disp.  Submerged  Disp.  +  Free  flood 
The  algorithms  chosen  for  the  weight  estimating 
module  take  standard  design  parameters  as  input,  where 
possible,  and  provide  an  output  in  terms  of  a  percentage 
of  Al  weights,  NSC,  or  submerged  displacement,  as 
appropriate.  These  percentages  are  then  combined  with 
the  discrete  weights  that  are  calculated  or  input  for 
selected  weight  groups  to  determine  the  Al  weight  total, 
lead,  NSC,  and  other  weight  parameters.  The  designer  has 
the  option  of  specifying  weights  for  groups  II,  IV,  VI, 
and  VII,  which  override  the  programmed  algorithms.  The 
following  section  discusses  the  chosen  algorithms  for 
calculating  each  weight  group.  Throughout  this  thesis, 
the  following  conventions  will  be  used  to  indicate 
mathmatical  operations: 

"+"  =  addition 

"-"  =  subtraction 

"*"  =  multiplication 

" / "  =  division 

ham  _  exponentiation 


GROUP  I : 

Weight  group  I  is  determined   as   a   percentage   of 
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NSC.  In  general, 

7.Wl=/(  depth,  material,  size). 

The  size  dependence  is  accomodated  by  calculating 
Wl  as  a  percentage  of  NSC.  The  algorithm  implemented  in 
the  program  is: 

7.W1  =  C1+C2*DEPTH. 

CI  and  C2  are  dependent  on  hull  material  used,  and 
are  determined  parametr i cal ly  from  data  on  past 
submarine  designs  C  4,5  11.  The  DEPTH  parameter  is  a  user 
input  of  maximum  operating  depth  in  feet.  The  accepted 
design  practice  of  designing  for  a  collapse  depth  of 
1507.  of  operating  depth  is  implemented  through  the 
choice  of  CI  and  C2. 

GROUP  I  I : 

Weight  group  II  is  calculated  directly,  in  tons, 
rather  than  as  a  percentage.  The  number  of  parameters 
involved  in  calculating  W2  is  greater  than  for  any  other 
weight  group. 

W2=/C horsepower ,  propulsion  type,  battery  type  &. 
capacity  ( capaci ty=endurance )  ] 

The  horsepower  variable  for  a  nuclear  plant  is 
simply  shaft  horsepower,  while,  for  a  non-nuclear  plant, 
two  separate  horsepowers  a.r&  needed;  shaft  horsepower 
(usually  the  horsepower  of  the  electric  final  drive 
motor)  and  charging  horsepower  (power  of  the  prime  mover 
used  to  drive  the  alternators  or  generators   to   provide 
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electrical  power  to  recharge  the  storage  batteries).  The 
weight  of  the  charging  prime  mover  is  a  strong  function 
of  technology  chosen;  le  diesel,  wankel,  etc. 

The  calculation  of  storage  battery  weight  for  a 
non-nuclear  plant  is  dependent  on  battery  type,  capacity 
in  kilowatt-hours  at  a  high  discharge  rate  ("sprint"), 
and  capacity  in  KW-H  at  a  low  discharge  rate  (endurance 
power).  The  options  provided  for  battery  type  in  this 
program  are  nickel-cadmium,  improved  lead-acid  (German 
VARTA  type),  and  "standard"  lead  acid  (U.S.  Trident 
type).  Constants  determined  from  the  energy  densities  of 
these  battery  types  at  the  two  discharge  rates  are  used 
to  determine  battery  weight  (the  larger  of  sprint  or 
endurance  rate  battery  weights  is  chosen  as  the  final 
battery  weight). 

Electric  motor  weight  is  determined  from  a 
multiplier  based  on  a  modern,  air-cooled  DC  motor.  The 
option  exists  to  provide  other  multipliers  based  on 
liquid  cooling,  superconducting  technology,  etc. 

For  the  non-nuclear  plant,  an  additional  option 
exists  to  provide  a  lumped  adjustment  to  weight  group  II 
to  provide  for  unconventional  propulsion  technology 
(fuel  cells,  Stirling,  low-power  nuclear,  etc.).  This 
final  adjustment  allows  the  designer  maximum  flexibility 
to  tailor  group  II  without  making  modifications  to  the 
program  code. 

For  a  "straight"  nuclear  plant, 
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W2=( C1*SHP) ./[ log (SUP)  DAC2 
where  SHP=shaft  horsepower 

CI  is  from  parametric  data. 
For  a  non-nuclear  plant, 
W2=WE+WM+WB+dW2 
where  WE=charging  prime  mover  weight 

WM=propulsion  motor  weight 

WB=battery  weight 

dW2=adJ ust men t  to  W2  (user  input) 
WE=C3*CHP 
where  CHP=charging  horsepower 

C3  is  dependent  on  prime  mover  type. 
WM=C4*SHP. 

WB=MaxC ( C5*KWHS ) , ( C6*KWHE ) ] 
where  KWHS= ca pa ci ty  for  "sprint" 

KWHE=capaci ty  for  endurance 

C5,C6  are    dependent  on  battery  type. 

GROUP  III: 

Although  weight  group  III  would  intuitively  appear 
to  be  a  direct  function  of  installed  generating 
capacity,  a  study  of  historical  design  data  C 5  1  reveals 
that  the  group  III  weights  are-,  in  fact,  very  closely 
approximated  by  a  straight  percentage  of  Al  weights. 
This  anomaly  is  best  explained  by  the  large  amount  of 
this  weight  group  attributable  to  the  electrical 
distribution  system   spread   throughout   the   submarine. 


-20- 


This  power  distribution  system  weight  is  a  direct 
function  of  the  size  of  the  submarine,  since  the  amounts 
of  cable,  connectors,  and  distribution  panels  vary  with 
size  of  the  vessel.  Weight  group  III,  consequently,  is 
calculated  as  a  percentage  of  Al  weights: 

y.W3=ci. 

GROUP  IV: 

Determination  of  group  IV  weights  is  v^r\j  difficult 
to  implement  in  an  "automatic"  mode.  The  variety  of 
combat  system  components  and  great  variation  in 
equipment  weights  makes  the  selection  of  an  algorithm 
extremely  difficult,  as  most  schemes  provide  accurate 
results  for  a  limited  number  of  combat  suites,  at  best. 
It  is  therefore  highly  recommended  that  the  designer 
specify  a  group  IV  weight  based  on  his  off-line 
determination  of  the  desired  combat  systems  and 
associated  weights.  To  permit  a  rough  approximation  is 
acceptable  for  an  early  design  iteration,  the  program 
contains  an  algorithm  based  on  historical  data  as  a 
percentage  of  Al  weights  for  either  fast  attack  or 
ballistic  missile  submarines: 

'/.W4=Ci   i  =  l,4 

where  Ci  is  based  on   submarine   type   (SS,   SSE, 
SSN,  or  SSEN). 


GROUP  V: 

Weight  group  V  is  accurately  determined  as  a 
function  of  submarine  size,  since  the  auxiliary  systems 
required  for  a  vessel  are  directly  proportional  to 
displacement  and  internal  volume.  Group  V  is  thus 
calculated  as  a  percentage  of  Al  weights: 

7.W5  =  Ci   i  =  l,2 

where  Ci  is  based  on  submarine   type   (attack   or 
ballistic  missile). 

GROUP  VI : 

Group  VI  weights  are  a  function  of  crew  size,  and, 
to  some  extent,  typical  mission  duration.  The  multiplier 
for  group  VI  weight  determination  is  based  on 
habitability  standards  commensurate  with  the  last  two 
classes  of  attack  and  ballistic  missile  submarines, 
respectively.  The  weights  a.re  directly  calculated  in 
tons  : 

W6=C1*NP 

where  NP=number  of  personnel  in  the  crew. 


GROUP  VII : 

The  estimation  of  armament  weights  is  nearly  as 
difficult  as  the  task  of  determining  group  IV  weights, 
since  the  possible  combinations  of  existing  and  future 
weapons  and  associated  launchers  are  infinite.  Again,  it 
is  highly  recommended  that  the  designer  specify  a   value 


for  weight  group  VII  based  on  his  off-line  determination 
of  the  weapons  payload.  As  was  provided  in  the  case  of 
combat  systems,  a  rough  estimate  of  group  VII  is 
available  from  a  percentage  of  Al  weights.  The 
percentage  constant  wa^s  determined  from  historical  data, 
and  separate  values  are  provided  for  attack  and 
ballistic  missile  submarines. 

y.w7=ci . 

OTHER  WEIGHTS: 

The   amounts   of   lead,   variable    load,    reserve 
buoyancy,  and  free  flood  are    provided   as   fractions   by 
the  designer.  The  following  ranges  are    appropriate: 
Lead  (.07-.ll)*Al 

Variable  Load       ( . 04-. 07 ) *NSC  (nuclear) 

( .08-. 19)*NSC  (non-nuclear ) 
Reserve  Buoyancy    ( . 12-. 15 ) *NSC 

Free  Flood  (.  05-.  1  )*(  Submerged  Displ..  ) 

The  operator,  however,  may  specify  his  own  values 
as  desired. 

CALCULATION  OF  WEIGHT  SUMMARY 

When  either  percentages  or  actual  values  have  been 
calculated  or  specified  for  each  of  the  weight  groups, 
the  calculations  for  the  weight  summary  are  processed. 
At  a  minimum,  actual  weight  values   in   tons   have   been 
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deter mined  for  weight  groups  II  and  VI.  Additionally, 
the  designer  may  have  specified  weights  for  weight 
groups  IV  and  VII.  As  an  example  of  the  algorithms  used 
for  final  weight  calculations,  assume  that  only  W2  and 
W6  have  been  determined  in  tons.  The  other  permutations 
of  final  weight  calculations  are  similar,  taut  allow  for 
dicrete  weights,  vice  percentages,  for  either  W4,  W7,  or 
both.  For  the  case  where  only  W2  and  W6  are  known,  the 
sequence  is  as  follows: 

(All  "'/."     figures  used  are    converted  to  fractions; 
7.X=7.X/100) 

K 1  =  (  7.W  1  +  7.W  1  *7.LD  )  /  (  1  -7.VL  )  +  7.W3+7.W5 

where  7.LD=lead  fraction 

7.VL=var  iable  load  fraction. 

A  1  =  (  W2+W6  )  /  (  1  -K  1  -7.W4-7.W7  ) 

W3=7.W3*A1 

W4=7.W4*A1 

W5=7.w5*Al 

W7=7.W7*A1 

Lead=7.LD*Al 

A=A1+I_ead 

NSC=A/(  1-7.VL) 

Variable    Load=7.VL*NSC 

Main    Bal  last=7.RBY*NSC 

where  7.REY=f  ract  ion  of  reserve  buoyancy 

Submerged  Disp. =NSC+Main  Ballast 

Free  Flood=7.FF*Submerged  Disp. 


where  7.FF=free  flood  fraction 

Envelope  Disp. =Submerged  Disp.+Free  Flood. 

The  calculated  weight  summary  is  displayed  on  the 
screen  with  a  summary  of  user  inputs,  at  which  time  the 
designer  may  accept  or  reject  the  results.  If  the 
results  are  rejected,  the  program  loops  back  to  the 
input  section  for  the  weight  module  and  another 
iteration  is  started.  If  the  results  ar&  accepted,  they 
are  written  into  an  output  file  for  record  purposes,  and 
various  input  and  calculated  parameters  ar&  written  into 
several  intermediate  files  for  use  in  later  program 
modules  as  the  design  is  completed. 

Because  of  the  high  degree  of  importance  attached 
to  the  calculated  weight  data,  future  refinement  of 
program  algorithms  should  include  a  priority  effort  to 
accomplish  two  major  tasks: 

(1)  Better  accuracy  of  calculated  weights. 

(2)  Further   subdivision   of   weights   (at   least 
partially  to  the  three-digit  SWBS  level). 
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PRINCIPAL  CHARACTERISTICS 
AND  ENVELOPE  GEOMETRY 

PRINCIPAL  CHARACTERISTICS 

The  results  from  the  weight  module  provide  the 
first  important  input  to  the  determination  of  the 
remainder  of  principal  characteristics,  which  ar& 
geometric  in  nature.  Since  this  thesis  implements  a 
weight-driven  submarine  design  process,  the  next  design 
task  is  to  create  an  envelope  of  the  appropriate  volume 
to  contain  the  envelope  displacement: 

Envelope  Volume=35*Envelope  Displacement. 
The  product  envelope  defined  by  this  module  of   the 
program  will  exhibit  the  following   geometric   principal 
characteristics    (abbreviations     and     symbols     in 
parentheses  a.re    program  variable  names): 

Length  Overall  (LOA) 

Length  of  Entrance  (Lf ) 

Length  of  Run  (La) 

Mid  Body  Length  ( Lm ) 

Diameter  (D) 
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Length /Diameter  Ratio  (L/D) 
Prismatic  Coefficient  (Cp) 
Eare  Hull  Surface  Area  ( SF ) . 
In  addition,  once  envelope  geometry  is  defined,  the 
following  calculated  hydrostatic  characteristics  may   be 
found : 

Submerged  Longitudinal  Center  of  Buoyancy 

(LCESUB) 
Draft  a  Normal  Surface  Displacement  (TNSC) 
Longitudinal  Center  of  Euiyancy  a  NSC  (LCBNSC). 
These  parameters,  then,  when  coupled  with  the  input 
parameters  delineated  for  the  weight  estimation   module, 
constitute  a  set  of  principal   characteristics   for   the 
design . 

DEFINING  THE  ENVELOPE 


A  number  of  algorithms  have  been  suggested  to 
calculate  the  geometry  of  the  envelope  for  submarines. 
These  range  from  a  " cut-and-paste"  approach  using  past 
designs    [2]    to    ver^  sophisticated     polynomial 

determination  of  offsets  C6D.  In  between  these  extremes 
Are  a  number  of  methods  utilising  simple  thumbrules  and 
parametric  relationships  among  Cp,  D,  L/D,  and 
displacement  [4,7]. 

Both  the  cut-and-paste  and  thumbrule  approaches  are 
somewhat  imprecise  in  allowing   the   designer   to   fully 
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control  the  geometry  of  the  resulting  hull  form.  The 
polynomial  methods,  while  very  precisely  controlled, 
require  the  designer  to  specify  design  parameters  which 
ar&  not  couched  in,  or  clearly  related  to,  traditional 
naval  architecture  vocabulary.  The  most  desirable 
algorithm,  then,  will  allow  strict  control  of  the 
resulting  geometry  while  accepting  input  parameters 
which  ar&    common  naval  architectural  terms. 

CHOSEN  ENVELOPE  GEOMETRY  ALGORITHM 

In  choosing  the  algorithm  implemented  in  this 
thesis,  two  restrictions  were  accepted  for  the  sake  of 
simplicity.  The  envelope  shape  is  confined  to  circular 
cross-sections,  and  an  analytically  smooth  envelope 
outline  is  calculated.  Figure  4.1  demonstrates  the 
nomenclature  and  subdivision  of  the  envelope  calculated 
with  the  chosen  algorithm. 

The  designer  is  prompted  for  the  following 
geometric  input  parameters: 

Length  of  Entrance 

Length  of  Run 

Forebody  Prismatic  Coefficient  (CpF) 

Aft  Body  Prismatic  Coefficient  (CpA) 

Maximum  Diameter  (D). 
The  algorithm  then  defines  an   envelope   consisting 
of  a  hemi-el 1 i pt ical  forebody,   a   cylindrical   midbody, 
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and  a  parabolic  afterbody.  The   chosen   algorithm   is   a 
modification  of  one  suggested  by  Jackson  [4]. 


Figure  4.1:  Nomenclature  for  determination   of   envelope 
geometry. 
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The  coefficients  Nf  and  Na  are  essentially  fullness 
coefficients  for  the  fore  and  aft  bodies,  respectively. 
In  order  to  facilitate  input  in  traditional  naval 
architectural  terms,  these  fullness  coefficients  are 
calculated  from  the  values  specified  by  the  designer  for 
prismatic  coefficients  of  the  fore  and  aft  body  sections 
of  the  envelope: 

Nf=  71.0477*CpFA5  +65 . 81 07*CpFA4  -483. 3037*CpFA3 
+587. 4137*CpFA2  -281 . 7224*CpF  +49.5876 

Where  Nf  is  the  forward  fullness  coefficient. 

Na=    379.6546*CpAA5  -938 . 4708*CpAA4 
+944.8853*CpAa3  -47 1 . 0B72*CpAA2  +1 19. 1465*CpA  -11.3454 

Where  Na  is  the  aft  fullness  coefficient. 
These  polynomial  relationships  are  valid  for  values 
of  CpF  from  .5383-. 8720,  and  for  CpA  from  .3333-. 7111. 
These  ranges  substantially  exceed  those  for  actual 
submarine  designs.  The  relationships  were  derived  by 
using  known  (hand-verified)  equivalencies  of  prismatic 
and  fullness  coefficients  to  set  up  matrix  equations 
with  the  polynomial  coefficients  as  the  unknown  column 
vector.  Gaussian  elimination  was  then  used  to  solve  for 
the  unk own  coefficients.  Additional  hand-calculated 
equivalencies  were  then  utilized  to  verify  that,  within 
the  specified  ranges,  the  polynomial  relationships  yield 
results  correct  to  the  fourth  decimal  place. 

The  "X"  coordinate  for  the  overall  envelope  runs 
from  0  at  the  forward  perpendicular  to  LOA   at   the   aft 
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perpendicular .   In   the    formula    for    the    forebody 
hem i -el  1  ipse,  a  local  X  coordinate  <Xf)  runs  from   0   at 
the  end  of  the  entrance  length  (Lf)  to  Lf  at  the  forward 
per pen di  cular : 
Xf  =  Lf-X 
In  the   aft   body   parabolic   section,   a   local   X 
coordinate  is  again  used: 
Xa=  X-(Lf+Lm) 
Note  that  Lm,  the  mid  body  length,  is  a   calculated 
parameter,  and  not  an  input.  The  calculation  of  Lm   will 
be  described  later  in  this  chapter. 

The  hull  radius,  r(X),  for  the  elliptical   forebody 
is  calculated  as  follows: 

r(X)=  R*[ ( l-(Xf /Lf ) ANf >A( 1/Nf ) ] 
where  R=  D/2 
The  hull  radius  within  the  mid  body  is  simply: 

r (X )=  R 
The  hull   radius   within   the   aft   body   parabolic 
section  is: 

r(X)=  R*C l-(Xa/La) ANaD 
The  offset  r(X)  is  calculated  for  each  of  21 
stations  evenly  spaced  along  the  hull  length.  Obviously, 
r(X)  at  both  the  forward  and  aft  perpendiculars  is  zero. 
When  the  offsets  for  the  21  stations  have  been 
determined,  the  sectional  ^rea  of  each  station  is 
calculated : 

SA(X)=  Pi*r(X)A2 
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The  station  sectional  areas  are  then  integrated 
along  the  length  to  determine  the  volume  enclosed  by  the 
envelope : 


=  /' 


'LOA 

/  SA' 
y0 


This  and  all  other  integrations  in  this  thesis  are 
done  numerically,  using  Simpson's  Rule. 

On  the  first  pass  through  this  process,  the  mid 
body  length  is  zero,  and  LOA=  Lf+La.  The  resulting 
volume,  then,  is  the  total  volume  for  the  fore  and  aft 
body  sections  alone.  To  determine  the  proper  mid  body 
length,  the  volume  determined  above  is  subtracted  from 
the  required  envelope  volume,  and  Lm  is  calculated  as 
the  length  of  a  cylinder  of  radius  R  containing  the 
remaining  volume: 

DV=  Envelope  Volume— VOL 
Lm=  DV/(Pi*RA2) 

If  DV  is  negative,  an  error  message  is  printed,  and 
the  designer  is  prompted  for  new  geometric  input  values. 

After  Lm  is  calculated,  the  final  value  of  LOA  is 
determined : 

LOA=  Lf+Lm+La 

The  sectional  areas  for  the  resulting  hull  form  are 
calculated  at  the  new  21  stations,  after  which 
circumferences  and  bare  hull  surface  area  are 
determined : 

Circ(X)=  2*Pi*r(>:) 
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•LOA 

Circ(X)  dX 


C 


The  overall  prismatic  coefficient  is  calculated   by 
comparing  the  envelope   volume   with   the   volume   of   a 
cylinder  (Vcyl)  of  length  LOA: 
Vcyl=  Pi*RA2*L0A 
Cp=  Envelope  Volume/Vcy 1 
These  calculated  parameters,   when   taken   together 
with   the   input   values,    constitute    the    geometric 
definition  of  the  envelope. 

CALCULATION  OF  EA5IC  HYDROSTATICS 

In  order  to  later  balance  the  submarine,  it  is 
important  to  accurately  determine  the  longitudinal 
centers  of  buoyancy  in  the  submerged  (LCESUB)  and  normal 
surface  (LCENSC)  conditions.  Calculation  of  LCENSC 
requires  the  determination  of  the  draft  at  NSC.  LCE7s  in 
both  conditions  are  calculated  by  integrating  the  first 
longitudinal  moment  of  the  sectional  areas. 

Using  the  sectional  areas  for  the  fully-submerged 
hull  as  calculated  above,  LCESUE  is  determined: 

rLOA 
LCBSUE=   /  X*SA(X>  dX  /VOLENV 


There    are  several    possible    approaches     to 

determining  the  draft  at  NSC.  When  the  calculations  are 
done  manually,  it  is  customary  to  choose  at  least   three 
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points,  commonly  at  drafts  equal  to  D,  D/2,  and  D+R/2, 
and  determine  the  volumes  and  resulting  displacements 
for  each  of  these  drafts.  A  curve  of  displacement  versus 
draft  is  then  constructed,  after  which  the  draft 
corresponding  to  the  NSC  displacement  is  read  from  the 
curve.  While  this  method  is  expeditious  and  acceptable 
for  hand  calculations,  relatively  large  errors  may 
result  when  the  designer  is  not  experienced  enough  to 
draw  the  proper  curve  through  the  three  points.  Use  of 
additional  points  improves  the  accuracy,  but  the  manual 
calculations  are    tedious  and  time-consuming. 

With  the  speed  and  accuracy  of  the  computer 
available,  a  much  more  precise  determination  of  draft  at 
NSC  (TNSC)  is  possible.  A  potentially  important  factor 
that  is  often  overlooked  in  manual  methods  is  the  volume 
of  the  free  flood,  which  is  contained  in  the  envelope, 
and  may  be  a  signifigant  portion  (2-10'/.)  of  the  envelope 
volume.  Study  of  existing  designs  indicates  that 
approximately  807.  of  the  total  free  flood  volume  is 
submerged  at  TNSC.  Consequently,  the  "target" 
displacement  in  determining  TNSC  is  not  simply  NSC,  but 
rather  includes  this  portion  of  the  free  flood: 
NSCVOL=  35*(NSC+.B*Free  Flood) 

(for  807.  free  flood).  Where 

NSCVOL  is  the  desired  volume  at  TNSC. 

The  method  chosen   to   determine   NSC   draft   is   a 
bracket   and   halve   iteration   on   the   variable   TNSC, 
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beginning  at  TNSC=  3*R/2.  The  allowable  error  in  NSC 
volume  is  specified  as  0.5/1  in  the  current  program.  This 
tolerance  may  be  decreased,  with  additional  iterations 
required  for  convergence.  The  value  of  0.5/1  was  chosen 
as  sufficient  for  the  conceptual  design  phase  while 
minimizing  the  number  of  iterations  required.  The 
algorithm  employed  is  as  follows: 

DV=  . 005*NSCVOL 

NPLUS=  NSCVOL+DV 

NMINUS=  NSCVOL-DV 

DR=  .5*R 

TNSC=  R+DR 
For  each  trial  TNSC,   the   value   of   each   station 
radius,  r(X),  is  compared  to  DR: 

If  (DR)>r(X),  then  r(X)  and  SA(X)  remain 
unchanged. 

If  (DRXr(X),  then: 

SA(X)=  SA(X)  -Cr(X)A2  *ARCOS ( DR/r ( X ) )  -DR 
*(r <  X) A2  -DRA2) A.5  3 

(Zero  trim  is  assumed). 

rlOA 

VOL=   /  SA(X)  dX 


If  VOL>NMINUS  and  VOL<NPLUS,  then  the  correct  TNSC 
is  determined.  If  VOL<NMINUS,  then  DR  is  increased  by 
.5*DR.  Similarly,  if  VOL>NPLUS,  then  DR  is  decreased  by 
the  same  increment.  In  either  case,  the  algorithm  is 
repeated  until  the  value  of  VOL  is  acceptable.  When   the 
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correct  TNSC  15  found,  LCENSC  ic  calculated: 

rlOA 

LCBNSC=   /  X*SA(X)  dX  /NSCVOL 
yB 

The  calculated  results  are  displayed  on  the  screen 
for  the  designer's  approval.  If  the  results  are  deemed 
unsatisfactory,  the  designer  may  reject  the  envelope 
geometry  and  return  to  the  input  section  for  another 
envelope  design  iteration. 

If  the  results  are  accepted,  all  of  the  specified 
input  values  and  calculated  data  are  written  to  an 
output  file  for  record  purposes,  and  appropriate  data 
for  other  program  modules  are  written  to  passing  files 
for  later  use. 

PLOTTING  THE  ENVELOPE 

An  accurate  plot  of  the  envelope  is  imperative, 
since  the  outline  displayed  on  the  monitor  will  guide 
the  designer  in  his  interactive  design  of  the  pressure 
hull.  The  21  stations  determined  thus  far  are 
insufficient  to  fit  an  accurate  curve,  since  the 
arbitrary  equally-spaced  stations  may  not  fall  at  the 
points  of  greatest  hull  curvature,  or  at  the  junctions 
of  the  mid  body  with  the  fore  and  aft  body  sections. 

To  adequately  portray  hull  curvature  and 
transitions  at  section  Junctions,  the  program  calculates 
132  unequally-spaced  offsets  to  be  used  in  plotting   the 
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envelope.  Offsets  are  most  closely  spaced  near  the 
forward  perpendicular  and  at  the  transition  points  of 
the  hull.  These  points  are  then  used  to  fit  a  Bezier 
spline  curve  which  forms  the  profile  outline  of  the 
envelope.  The  plot  is  automatically  scaled  to  make  full 
use  of  available  screen  area-,  and  both  vertical  and 
horizontal  scales  are  drawn  as  a  reference  for  the 
designer . 

The  two-dimensional  envelope  outline  is  rotated  to 
form  a  surface  of  revolution,  and  a  mesh  pattern  is 
superimposed  on  the  surface  to  aid  visualization  of  the 
resulting  three-dimensional  hull  form.  In  addition,  an 
isometric  view  is  displayed  below  the  profile  for 
perspect  i ve. 

The  resulting  display  of  two  views  of  the  outer 
hull  constitutes  the  graphic  output  of  the  envelope 
geometry  module.  The  designer  is  now  ready  to  proceed  to 
the  speed  and  power  calculations. 
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SPEED  AND  POWER 


GENERAL  DISCUSSION 

The  calculation  of  resistance  for  the  envelope 
geometry  developed  in  the  previous  module  is  a 
straigh tf oruard  analytic  process.  The  formulae  employed 
Are  the  familiar  drag  estimation  relationships  from 
basic  hydrodynamics.  A  proposal  for  future 
implementation  would  be  the  addition  of  a  propeller 
design  module  to  accurately  predict  the  propulsive 
coefficient,  in  conjunction  with  the  envelope  geometry., 

To  facilitate  endurance  calculations  for 
non-nuclear  designs,  the  calculation  of  total  electrical 
load  at  each  speed  is  included.  This  total  load  figure 
must  include  all  non-propulsion  loads  (hotel  and  combat 
system)  as  well  as  those  required  to  propel  the 
submar  ine. 

With  the  data  passed  from  the  envelope  geometry 
module,  an  estimate  of  appendage  drag  may  be  made.  This, 
when  combined  with  hull   form   and   surface   area   data, 
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completes  the  input  to  the  speed  and  power  algorithm. 
ESTIMATION  OF  APPENDAGE  DRAG 

Study  of  existing  submarines  reveals  that  appendage 
drag  comprises  thirty  percent  or  more  of  the  total  drag 
on  the  submarine  for  most  designs.  In  the  later  stages 
of  the  design  process,  detailed  consideration  should  be 
given  to  the  configuration  and  size  of  each  appendage  to 
improve  the  accuracy  of  resistance  estimates.  At  the 
stage  of  conceptual  design,  however,  details  concerning 
appendages  are  often  uncertain.  Many  designs  are 
formulated  without  final  selection  of  control  surfaces, 
for  example.  The  required  shape,  location  and  sizes  of 
surfaces  for  cruciform,  X,  and  inverted  Y  sterns  may 
exhibit  some  variance,  and  the  final  configuration  may 
not  be  chosen  until  much  later  in  the  design  process. 

In  view  of  this  uncertainty,  a  reliable  rough 
estimate  of  total  appendage  drag  is  desirable  for  the 
conceptual  design.  This  thesis  incorporates  the  method 
suggested  by  Bukalov  LQ1.  Eukalov  studied  numerous 
actual  submarine  designs,  and  developed  an  algorithm 
based  on  this  parametric  study.  This  algorithm  may  be 
reduced  to  a  single  formula  dependent  on  length  overall 
(LOA)  and  maximum  diameter  (D).  The  appendage  drag  (Da) 
is  then  calculated  as  follows: 

Da=  1 .09065E-3*LOA*D+1 1.25 
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CH05EN  ALGORITHM 

The  basic  relationships  implemented  in  this  program 
module  may  be  found  in  numerous  introductory  texts 
[9,10].  The  interactive  input  from  the  terminal  consists 
of  the  following  items: 

Propulsive  Coefficient  (PC) 
Non-propulsion  Loads  in  KW  (HL) 
The  remainder  of  the  required  input  is  passed   from 
the  envelope  geometry  module,  and  consists  of: 
Prismatic  Coefficient  (Cp) 
Length  Overall  (LOA) 
Maximum  Diameter  (D) 
Eare  Hull  Surface  Area  ( SF ) 
The  speed  and  power  data  is  calculated   and   stored 
in  an  array  indexed  on  speed  in  knots  (Vk)  from  zero   to 
forty  knots.  After  appendage  drag  is   estimated   in   the 
manner  previously  described,  the  calculations  proceed  in 
an  iterative  loop: 

Re ( Vk ) =  ( Vk* 1 . 689*L0A ) / 1 . 27908E-5 
Where  Re  is  the  Reynolds  Number. 
Cf=  .075/C log(Re)-23A2 

Where  Cf  is  the  frictional  drag  coefficient, 
calculated  by  the  ITTC  convention. 

0=  Cf  *C  1  .  5*  (  D/LOA  )  A  1  .  5+7*  (  D/LOA  )  A3  ] 
+. 002* ( Cp-. 6) 
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Where  Cr  is  the  residual,  or  form,  drag 
coefficient,  calculated  after  Jackson  CAD,  and  the 
last  term  is  a  correction  for  submarines  with 
considerable  parallel  mid  body. 
Ct=  Cf+Cr+. 00025 

Where  Ct  is  the  total  drag  coefficient  for  the 
bare  hull,  and  the  final  term  is  the  correction 
factor  from  tow  tank  studies. 

HP( Vk)=  ( l/PC)*.00S72*VkA3*( Ct*SF+Da) 
Where  HP(Vk)  is  the  drag  in  horsepower  for  the 
given  speed,  Vk. 

KW(Vk)=  .  7A57*HP(v,k)+HL 

Where  KW(Vk)  is  the  total  load  in  kilowatts, 
including  n on -pro pulsion  loads,  for  the  speed,  Vk. 
These  calculated  results  for  drag  in  horsepower  and 
total  load  in  kilowatts  ar&  displayed  on  the  terminal 
screen  for  designer  inspection.  From  these  results,  the 
designer  may  determine  the  maximum  speed  of  the 
submarine  for  the  installed  shaft  horsepower,  as  well  as 
the  submerged  endurance  at  any  speed  for  installed 
energy  storage  capacity.  If  desired,  the  module  may  be 
re-run  for  various  values  of  propulsive  coefficient  and 
non-propulsion  loads,  until  a  satisfactory  result  is 
obtained. 

When  the  designer  accepts  the  calculated  values, 
the  results  are  written  into  an  output  file  for 
hard-copy  record  purposes. 
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PRELIMINARY  BALANCE 


GENERAL  DISCUSSION 

When  the  designer  has  completed  the  speed  and  power 
module,  he  will  proceed  to  the  design  of  the  pressure 
hull,  which  is  a  process  covered  separately  in  a 
concurrent  Ocean  Engineer  thesis  by  Marvin  Meade  [11  D. 
With  a  graphic  display  of  the  envelope  and  pressure  hull 
to  serve  as  a  visual  aid,  the  designer  will  enter  the 
subject  module  for  this  chapter  to  accomplish  a 
longitudinal  weight  and  buoyancy  balance. 

The  purpose  of  the  preliminary  balance  is  to 
determine  the  locations  of  lead  and  main  ballast  tanks, 
in  order  to  confirm  the  feasibility  of  the  pressure  hull 
design.  The  required  lead  (Pb)  longitudinal  center  of 
gravity  (LCG)  is  achieved  by  balancing  the  longitudinal 
moments  of  all  weights  except  the  main  ballast  tanks  to 
create  an  attitude  of  zero  trim  in  the  normal  surface 
condition  (NSC). 

Following  the  determination  of  lead   placement,   an 
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additional  moment  balance  is  computed  to  locate  the  main 
ballast  tank  LCG  in  order  to  achieve  zero  trim  in  the 
submerged  condition. 

In  addition  to  the  longitudinal  balance,  a  vertical 
weight  balance  is  computed  to  provide  an  estimate  of 
transverse  stability  while  submerged.  This  estimate  is 
in  the  form  of  EG,  the  distance  from  the  vertical  center 
of  buoyancy  to  the  vertical  center  of  gravity  of  the 
submarine.  This  parameter  should  have  a  value  of  about 
one  foot  for  an  optimum  design.  To  simplify 
calculations,  no  trim  is  allowed  in  the  surfaced 
condi  t  ion . 

CALCULATION  PROCESS 

The  first  set  of  input  data  required  for  balancing 
the  submarine  is  passed  from  previously  completed 
modules.  The  data  which  is  read  in  from  earlier 
calculations  is  as  follows: 

All  16  weight  categories  from  the  Weight  Module 
(W1-W16) 

Longitudinal  Centers  of  Buoyancy  at  NSC  and 
submerged  (LCBNSC  and  LCBSUB) 
Maximum  Envelope  Radius  (R) 
Length  Overall  ( LOA ) . 
The  first  set  of  keyboard  input  is  the  LCG  and   VCG 
information  for  the  seven  primary  weight  groups  and   the 
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variable  load: 

LCG< 1 )-LCG(7) 

VCG< 1 )-VCG<7) 

LCG(VL) 

VCG(VL) 
Although  detailed  arrangements  have  not,  as  yet, 
been  addressed  in  the  design  process,  the  designer  must 
have  given  enough  thought  to  rough  arrangements  to 
provide  these  location  data  for  the  aggregate  weight 
subdivisions.  If,  when  detailed  arrangements  are 
completed  at  some  later  time,  weight  group  LCG  and  VCG 
locations  are  signifigantly  different  from  the  input 
provided  for  this  module,  the  balance  process  must  be 
repeated. 

The  program  now  calculates  arrays   of   longitudinal 
and  vertical  moment  data: 

LM< l )=  LCG( l )*W< i ) 

VM( i )=  VCG( i )*W(i ) 

Where  LM  and  VM  are    the  longitudinal  and  vertical 
moments  of  the  weight  groups,  respectively. 


LCG(Cond.  A-l)=  >  LM(i)  /  /    W(i) 


VCG ( Con d.  A-l)=  >  VM ( i )  /     /    W(i) 


The  designer  is  now  given  a  screen  display   of   the 
total  weight  of  lead,  and  is  asked  to  input  the   desired 
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amount  of  that  total  to  be  used  as  margin  lead  (ML). 
This  lead  will  arbitrarily  be  located  at  LOA/2 
longitudinally  and  D/2  vertically.  The  remainder  of  the 
lead  is  considered  to  be  trim  lead  ( TL) .  The  required 
LCG  for  trim  lead  is  then  calculated: 

LCG(TL)=  CW(NSC)*LCBNSC  -LM(Cond.  A-l )  -LM ( VL) 

-ML*L0A/2  D/SL 

Similarly,  the  actual  VCG(Lead): 

VCG(Lead)=  (TL*6  +ML*R ) /W ( Lead ) 

Where  6  feet  above  the  keel  is  arbitrarily  chosen 

for  VCG(TL)  to  make  placement  in  the  circular  cross 

section  of  the  ballast  tanks  feasible. 

The  designer  is  then  asked  whether  the  required 
LCG(SL)  is  feasible.  If  the  response  is  affirmative,  the 
program  continues  on  to  main  ballast  tank  ( MBT ) 
calculations.  If  not,  the  designer  is  asked  for  new  LCG 
and  VCG  data  for  the  weight  groups  and  variable  load.  He 
is  prompted,  in  each  case,  with  the  previously  entered 
LCG  and  VCG  for  each  category.  Of  course,  if  reasonable 
adjustments  in  weight  group  locations  cannot  lead  to 
feasible  stability  lead  LCG  (feasible  being,  generally, 
between  the  physical  centers  of  the  main  ballast  tanks), 
a  re-appraisal  of  the  pressure  hull  design  and  rough 
arrangements  is  necessary. 

When  the  lead  has  been  satisfactorily  located, 
calculation  of  required  MBT  LCG  is  accomplished 
(VCG(MBT)  is  arbitrarily  fixed  at  R): 
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LCG(NSC)=    CLMCCond.     A-l  )    +LM(Lead) 

+LM< VL)  :/W(NSC) 

LCG ( MET ) =    C  W ( Submerged ) *LCBSUE 

-LCG ( NSC ) *W ( NSC ) 3/W ( NET ) 

The  designer  is  again  asked  to  decide  whether  the 
calculated  required  LCG  is  feasible,  with  similar 
results  in  program  flow  depending  on  the  answer-. 
Feasibility  of  MET  location  is  dependent  on  the  desired 
geometry  of  physical  locations  for  METs  as  envisioned  by 
the  naval  architect.  The  continued  graphic  display  of 
envelope  and  pressure  hull  geometry  is  useful  in  aiding 
this  decision  process. 

The  program  now  sums  all  vertical  moments  for 
weight  groups,  lead,  variable  load  and  main  ballast  to 
arrive  at  an  overall  VCG  in  the  submerged  condition 
(VCGSUE).  An  estimate  of  submerged  stability  may  be 
calculated  and  displayed: 
BG=  R-VCGSUB 

Note  that  the  vertical  center  of  buoyancy  is 
assumed  to  be  at  R  (D/2)  for  the  submarine,  which  is 
basically  a  body  of  revolution.  Any  error  due  to 
buoyancy  of  appendages  will  be  small  and  conservative, 
since  off-center  buoyancy  items  (sail,  etc.)  are 
generally  above  the  geometric  center  of  the  envelope. 
Again,  the  designer  is  asked  to  accept  or  reject  the 
calculated  stability,  and,  again,  rejection  results  in 
looping  back  to  the  input  section  of  the  program  module. 
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When  all  calculated  values  s,r&  deemed  acceptable, 
the  input  data  and  results  are  written  to  an  output  file 
for  record  purposes,  and  the  locations  of  each  item's 
center  of  gravity  are  plotted  on  the  pressure  hull 
graphic  display.  The  designer  is  now  ready  to  proceed  to 
the  calculation  of  the  equilibrium  polygon. 
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THE  EQUILIBRIUM  POLYGON 


GENERAL  DISCUSSION 

One  of  the  most  important  sets  of  calculations  in 
the  submarine  design  process  is  the  determination  of  the 
equilibrium  polygon.  Any  submarine  experiences  changes 
in  weight  distribution  during  the  course  of  operations, 
as  items  in  the  variable  load  are  consumed.  A  system  of 
tanks  must  be  sized  and  located  in  such  a  way  that  this 
change  in  weight  distribution  and  amount  may  be 
compensated  to  allow  the  submarine  to  remain  in  proper 
balance  and  trim. 

A  detailed  list  of  the  variable  load  items  to  be 
considered  may  be  found  later  in  this  chapter.  In 
aggregate,  the  items  compose  major  groups  in  the 
categories  of  provisions,  ammunition,  stores,  fuel,  and 
other  consumable  fluids.  For  a  nuclear  or  other 
non-f ossi 1-f uel  submarine,  the  boat  will  always  become 
lighter  as  items  are  consumed,  and  the  compensation 
consists   of   adding   appropriate   amounts   of   variable 
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ballast  at  the  proper  locations.  In  the  case  of  a  diesel 
or  other  fossil-fuel  submarine,  a  paradox  exists  for  the 
designer.  As  fuel  is  consumed  from  the  fuel  ballast 
tanks  external  to  the  pressure  hull,  these  tanks  are 
compensated  with  seawater,  since  they  are  "soft"  tanks, 
and  must  remain  full  of  fluid  to  maintain  a  zero 
differential  pressure  across  the  tank  structure.  The 
seawater  used  for  compensation  has  a  higher  density  than 
the  consumed  fuel  which  it  displaces.  Consequently,  the 
submarine  actually  becomes  heavier  as  fuel  is  consumed. 

The  tool  which  is  used  to  calculate  the  required 
locations  and  capacities  of  the  compensating  tanks  is 
called  the  equilibrium  polygon,  which  is  a  plot  of 
weight  versus  moment  (fore  and  aft)  of  the  weight. 
Figure  7.1  is  an  example  of  the  equilibrium  polygon  for 
a  submarine.  The  points  within  the  polygon  represent 
extreme  operating  conditions  for  the  submarine  in  terms 
of  the  weight  and  moment  required  to  balance  the 
conditions.  The  solid  lines  forming  the  "polygon"  are 
plots  of  the  weight  and  moment  added  by  progressive 
filling  and  emptying  of  the  major  compensating  tank 
groups.  These  groups,  which  may  be  comprised  of  several 
tanks  each  (the  aggregate  weight  and  moment  for  the 
group  is  plotted),  are  the  forward  trim,  auxiliary,  and 
after  trim  tank  groups.  If  all  of  the  operating  points 
for  the  submarine  are  contained  within  the  polygon,  the 
submarine  may  be  safely   compensated   for   all   expected 
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operat lonal  conditions.  If  not,  the  designer  must  change 
the  location(s),  weight(s),  or  both,  for  the  appropriate 
group(s),  in  order  to  expand  the  polygon  to  include  all 
the  conditions. 

The  conditions  to  be  calculated  for  polygon  points 
are  specified  in  the  NAVSHIPS  Technical  Manual  ( NSTM ) , 
Chapter  9290  C12D.  The  extreme  conditions  are  more 
severe  than  an  actual  submarine  would  be  expected  to 
encounter,  and  represent  a  "worst  case"  situation  for 
weight  and  moment  balance.  These  extreme  conditions  may 
be  thought  of  as  representing  the  following  operational 
si tuat  ions : 

HEAVY  #1:  Following  a  short,  fast  patrol,  during 
which  no  ordnance  is  expended,  but  all  fuel  oil  is 
consumed. 

HEAVY  #2:  Same  as  the  previous  condition,  except 
that  only  the  fuel  oil  in  the  fuel  ballast  tanks  is 
consumed. 

LIGHT  #1  &  LIGHT  #2:  Two  variations  following  a 
short  patrol  with  no  fuel  consumed,  but  with  all 
ordnance  expended. 

HEAVY  FORWARD  #1 :  Following  a  patrol  where  only  aft 
ordnance  is  expended,  and  all  fuel  oil  is  consumed 
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from  the  forward  fuel  ballast  tanks 

HEAVY  FORWARD  #2::  Same  as  previous  condition, 
except  that  fuel  from  normal  fuel  oil  tanks  is  also 
consumed. 

HEAVY  AFT:  Following  a  patrol  with  only  forward 
ordnance  expended,  and  with  fuel  oil  from  the  aft 
fuel  ballast  and  normal  fuel  oil  tanks  consumed. 

CONDITION  N:  The  "normal"  full-load  patrol 
condition  at  the  beginning  of  a  patrol.  For 
submarines  with  FETs,  the  FETs  B.re    specified  as 
fully  ballasted  with  seawater. 

C 0 ND I T 1 0 N  M  :  T  h  e  same  as  con  d  i  1 1  o  n  N ,  e  x  ce p t  t  h  a  t 
submarines  with  FETs  carry  a  full  fuel  load  in  the 
FETs  (for  such  submarines*  this  condition  will  be 
1 i  g  h  ter  t  h an  con  d  i  t  i  on  N  )  . 


For  the  above  extreme  conditions,  water  density  is 
also  specified  within  typical  ocean  ranges,  so  as  to 
accentuate  the  worst  case  nature  of  the  operating 
condi  t  ions. 
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THE  VARIABLE  LOAD 

The  items  that  comprise  the  category  referred  to  as 

the  variable  load  &r&  numerous  and  diverse  in  nature?., 
C 1  as si.fi  C3. t  i on  of  the  i  terns  i s  ag a  i n  a i  ded  by  t h e 
guidelines  and  categories  specified  in  the  NAVSHIPS 
Technical  Manual. 

The  variable  load  is  subdivided  into  two  large 
sub-groups;  fixed  items,  not  expected  to  var^ 
significantly  during  a  patrol,  and  the  truly  variable 
items,  routinely  expended  as  the  submarine  operates.  The 
"fixed"  portion  of  the  variable  load  consists  of  the 
following  items: 

CREW  AND  EFFECTS:  Actual  weight  of  the  submarine's 
complement  of  personnel,  clothing,  and  associated 
personal  effects. 

SLEMs  OR  COMPENSATING  WATER:  The  weight,  for  a 
fleet  ballistic  missile  (FBM)  submarine,  of  the 
submarine  launched  ballistic  missiles  (SLEMs)  or 
the  compensating  water  in  their  absence.  SLEM 
compensating  water  exactly  matches  the  weight  and 
moment  of  the  missiles  themselves,  so  that  this 
category  may,  indeed,  be  considered  a  fixed  item. 
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SANITARY  TANKS:  The  weight  of  tanks  and  contents 
associated  with  waste  disposal. 

OXYGEN  CANDLES:  The  weight  of  these  devices,  used 
for  the  generation  of  oxygen  under  emergency 
conditions.  These  are    generally  a  very  small  weight 
i  tern. 

LUEE  OIL  IN  SUMPS:  That  weight  of  lube  oil 
contained  in  the  oil  sumps  of  equipment. 

FIXED  CLEAN  FUEL  OIL:  That  weight  of  oil  carried  in 
non-compensated  tanks  which  are    maintained 
essentially  full  (i.e.,  shield  tanks  on  nuclear 
submar  ines ) . 

As  can  be  seen  from  Reference  [12D,  there  are  other 
items  specified  as  fixed  load,  but  they  are  minor  in 
nature,  and  beyond  the  level  of  detail  required  for 
conceptual  design.  For  example,  the  item  "depth  control 
tanks"  in  the  NSTM  is  not  included.  This  omission  will 
be  explained  in  the  discussion  of  fuel  load  subdivision 
in  the  following  section. 

The  items  considered  variable  in  nature  are  as 
f ol lows : 

PROVISIONS  AND  STORES:  The  consumable  foodstuffs 
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(provisions)  and  spare  parts,  paper  supplies,  etc. 
(stores)  that  the  submarine  crew  will  use  during  a 
patrol . 

REVITALIZATION  OXYGEN:  Oxygen  used  for  life  support 
(atmosphere  control)  during  lengthy  submerged 
periods.  On  newer  submarines  (particularly 
nuclear),  a  method  of  oxygen  generation  often 
obviates  the  need  for  oxygen  banks. 

TORPEDOS,  MISSILES,  AND  AMMUNITION:  This  category 
is  self-explanatory.  All  expendable  ordnance  other 
than  SLEMs  is  included. 

WRT  TANKS:  This  item  is  the  "water  round  torpedo" 
tank  capacity,  used  to  compensate  for  the  presence 
or  absence  of  torpedos  in  the  tubes. 

RESERVE  ELECTROLYTE:  The  electrolyte  carried  to 
replenish  the  submarine  storage  battery. 

FRESH  UATER:  For  this  thesis,  this  is  a  composite 
category,  consisting  of  the  total  of  potable,  feed, 
and  battery  water. 

RESERVE  LUEE  OIL:  Lube  oil  carried  to  replenish  the 
sump  lube  oil  during  the  patrol. 
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FUEL:  Propulsion  fuel  for  the  submarine.  See  the 
discussion  that  fol lows  for  the  sub-categories  of 
fuel  and  their  significance  to  the  polygon 
computat  ions. 

Fossil-fuel  submarines  present  a  variable  load 
parodoK  for  the  designer.  Fuel  which  is  carried  external 
to  the  pressure  hull  (the  majority  of  the  fuel  load  for 
most  "conventional"  submarines)  must  be  replaced  by 
seawater  as  it  is  burned,  since  these  external  "fuel 
ballast  tanks"  (FETs)  ar&  not  constructed  to  withstand  a 
high   pressure   differential   (i.e.,   they    &r&  "soft 

tanks").  The  seawater  which  is  used  to  compensate  these 
tanks  is  heavier  than  the  fuel  it  replaces.  Thus,  as  the 
submarine  consumes  fuel,  it  gets  heavier,  rather  than 
1 ighter . 

To  provide  a  margin  of  safety  at  the  base  of  the 
equilibrium  polygon,  some  means  must  be  provided  to 
remove  weight  as  fuel  is  burned.  In  earlier  designs,  a 
"safety  tank"  was  provided,  of  sufficient  capacity  that 
when  it  was  emptied,  enough  water  weight  was  pumped 
overboard  to  compensate  for  the  added  seawater  in  the 
FETs  as  fuel  was  consumed.  The  safety  tank,  then,  was 
full  at  the  beginning  of  the  patrol,  and  would  be  empty, 
or  nearly  so,  when  all  fuel  had  been  removed  from  the 
FETs.  For  purposes   of   this   thesis,   the   safety   tank 
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concept  was  not  used,  since  the  weight  and  volume  of  the 
tank  and  it's  contents  is  a  non-productive  load  item. 
However,  a  means  to  compensate  for  the  FET  weight 
addition  is  necessary. 

The  method  chosen  for  implementation  in  this  thesis 
is  the  provision  for  a  variable  fuel  oil  (VFO)  tank.  The 
VFO  tank  is  internal  to  the  pressure  hull,  and  does  not 
require  seawater  compensation  as  fuel  is  consumed.  In 
actual  practice,  this  tank  may  be  ballasted  and  treated 
as  a  dual-purpose  tank,  acting  as  an  additional  trim  or 
auxiliary  tank.  For  computational  simplicity  within  the 
program  module,  however,  the  VFO  tank  is  treated  as  an 
uncompensated  clean  fuel  oil  tank.  An  operational  mode 
of  employment  is  assumed  whereby  the  fuel  in  the  VFO 
tank  is  consumed  first  during  the  patrol,  before  the 
FETs  are  utilized.  Thus,  the  size  of  the  VFO  tank  must, 
at  a  minimum,  be  sufficient  to  compensate  for  the 
difference  in  fuel  and  seawater  densities  for  the  FET 
capacity.  This  dictates  a  minimum  capacity  of  23"/.  of  the 
total  fuel  load.  The  designer  is  asked  to  input  the 
percentage  of  total  fuel  to  be  allocated  to  the  VFO  tank 
(or  tank  group,  as  arrangements  may  dictate).  If 
sufficient  internal  volume  is  available,  it  may  be 
advantageous  to  use  a  figure  greater  than  23"/.,  in  order 
to  provide  an  adequate  margin  at  the  base  of  the 
polygon . 
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CALCULATING  THE  POLYGON 

The  program  module  execution  begins  by  reading  the 
weight  group  values,  submerged  LCE,  and  LCG  of  the  main 
ballast  tanks  from  data  files  generated  by  the  weight 
estimating  and  balance  modules.  The  total  load  to 
submerge  the  submarine  is  calculated: 

Load  to  Submerge=  Submerged  Displacement-  Condition 

A  Weights. 

The  total  variable  load  weight  is  displayed,  and 
the  designer  is  prompted,  item  by  item,  for  the  weight 
and  LCG  of  each  variable  load  item.  As  each  weight  and 
LCG  are  entered,  the  balance  of  the  weight  in  the 
variable  load  account  is  displayed.  Items  which  are 
commonly  distributed  in  two  sub-items,  fore  and  aft,  are 
input  as  two  separate  weights  and  LCGs.  Examples  of  such 
items  are  torpedos,  missiles,  and  fuel.  The  designer 
should  have  given  some  thought  prior  to  this  input 
sequence  to  the  subdivision  of  the  variable  load.  The 
envelope  and  pressure  hull  are  displayed  at  the  top  of 
the  screen  during  this  stage  of  the  module,  to 
facilitate  estimating  the  LCG  for  each  item. 

The  last  item  to  be  entered  is  fuel.  The  balance  of 
the  variable  load  account  at  this  point  is  assumed  to  be 
all  fuel  (obviously,  the  balance  should  be  zero  for  a 
non-fossil-fuel   submarine).   The   total   fuel   load   is 
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displayed,  and  the  designer  is  asked  to  enter  the 
fraction  to  be  allocated  to  the  VFO  tank.  He  is  then 
asked  to  provide  LCGs  for  the  forward  and  aft  FETs  and 
the  VFO  tank. 

With  the  input  of  all  items  completed,  the  program 
displays  a  variable  load  summary  of  all  item  weights  and 
LCGs.  If  the  summary  is  satisfactory,  the  program 
proceeds  with  calculations;  if  not,  the  designer  is 
returned  to  the  input  sequence,  where  the  prompt  for 
each  item  now  includes  a  display  of  its  current  value. 
If  this  value  is  to  remain  unchanged,  the  operator 
simply  hits  the  "RETURN"  key;  otherwise,  a  new  value  may 
be  specified.  This  iterative  input  procedure  is  repeated 
until  the  designer  is  satisfied  with  the  variable  load 
summary . 

When  the  input  of  all  variable  load  items  is 
complete  and  accepted  by  the  operator,  the  program 
calculates  the  water  to  balance,  and  the  associated 
moment  of  that  water,  for  each  of  the  nine  specified 
polygon  points.  Note  that  conditions  M  and  N  will  be 
identical  for  nuclear  submarines.  Condition  M  is  only 
defined  for  submarines  with  FETs. 

A  computation  of  one  of  the  equilibrium  conditions 
would  proceed  in  the  following  manner.  Prior  to 
computing  the  individual  conditions,  the  arms  and 
moments  for  each  of  the  variable  load  items  ^r& 
determined : 
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Item  Arm=  Submerged  LCE  -Item  LCG 

Item  Moment =  Item  Weight  *  I  tern  Arm 

In  addition,  factors  of  proportionality  are 
determined  for  light  and  heavy  seawater.  For  each 
condition,  the  sum  of  all  variable  load  items,  in  the 
fractional  amount  specified  by  the  NAVSHIPS  Technical 
Manual,  is  computed.  Similarly,  the  total  moment  of  the 
specified  items  is  determined.  The  required  water  to 
balance  and  moment  to  balance  may  now  be  calculated: 

Water  to  Ealance=  Load  to  Submerge  -MET  -Sum  of 

Item  Weights 

where  MET=  Main  Ballast  Weight 

Moment  to  Ealance=  -MET  Moment  -Sum  of  Item  Moments 

If  required,  the  water  and  moment  to  balance  are 
multiplied  by  the  proportionality  factor  for  light  or 
heavy  seawater.  Each  of  the  nine  conditions  is 
calculated  in  a  similar  fashion,  resulting  in  arrays  of 
weight  values  in  tons  and  moments  in  foot-tons  (with 
negative  moment  values  indicating  an  aft  moment). 

Several  "housekeeping"  chores  are  performed  at  this 
point  in  the  module.  The  display  of  envelope  and 
pressure  hull,  and  associated  data  base  attributes,  are 
filed  under  the  part  name  MIT. HULLOUT,  for  later  recall 
to  plot  desired  views  of  this  graphic  screen.  The 
maximum  values  of  water  and  moment  to  balance  are 
examined  to  determine  the  required  scale  for  the  polygon 
Plot.  A  new  part  is  activated  under  the  name  MIT.POLYOUT 
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for  graphically  displaying  the  polygon.  Depending  on  the 
scale  required,  one  of  several  pre-constructed  "form" 
drawings  is  activated,  providing  the  axes  and 
appropriate  scale  information  for  the  polygon  under 
con si  derat ion . 

The  program  now  inserts  and  labels  each  of  the 
equilibrium  condition  points  on  the  graphic  display. 
This  gives  the  designer  a  clear  display  of  the  limits 
required  for  his  polygon  design.  An  important 
consideration  at  this  point  is  whether  an  adequate 
margin  remains  at  the  base  of  the  polygon.  Generally,  if 
the  lowest  weight  value  is  less  than  15  tons  from  the 
abscissa,  the  designer  should  seriously  consider 
re-structuring  the  variable  load,  as  necessary,  to 
achieve  this  margin.  This  may  be  accomplished  by 
returning  to  the  input  section  of  the  module.  An  example 
of  such  a  change  would  be  to  place  more  of  the  fuel 
load,  volume  permitting,  in  the  VFO  tank,  for  a 
fossil-fuel  submarine. 

The  operator  is  now  asked  to  provide  the  LCGs  of 
the  Forward  Trim,  Auxiliary,  and  After  Trim  tank  groups. 
He  is  then  prompted  for  the  capacities,  in  tons,  for 
each  of  the  tank  groups.  The  program  determines  the  arms 
and  moments  of  each  group  in  the  same  fashion  previously 
delineated  for  the  variable  load  items. 

Lines  are  now  inserted  on  the  display  to  represent 
the  successive  filling  of  Forward  Trim,   Auxiliary,   and 
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After  Trim  tanks,  fol lowed  by  the  successive  draining  of 
the  tanks  in  the  same  order.  These  lines  form  the 
polygon,  and  all  equilibrium  points  must  fit  within  its 
boundaries  if  the  tank  locations  and  capacities  are 
adequate.  Obviously,  a  polygon  whose  boundaries 
substantially  exceed  the  extent  of  the  equilibrium 
points,  while  safe,  represents  a  waste  of  internal 
volume,  and  should  not  be  accepted. 

The  designer  now  has  the  option  to  accept  or  reject 
the  polygon,  as  plotted.  If  it  is  rejected,  the  polygon 
boundaries  are  erased,  and  the  designer  may  specify  new 
tank  group  capacities,  LCGs,  or  both,  after  which  a  new 
polygon  is  plotted.  This  process  may  be  repeated  as  many 
times  as  are    required  to  achieve  an  acceptable  polygon. 

When  the  polygon  results  are  satisfactory,  the 
polygon  is  plotted  on  the  electrostatic  plotter,  the 
part  MIT.POLYOUT  is  deleted,  and  the  variable  load 
summary,  tank  locations,  and  tank  capacities  are  written 
to  an  output  file  for  a  permanent  record  of  the  final 
input  data  set.  One  iteration  of  the  submarine 
conceptual  design  process,  as  supported  by  this  thesis, 
is  complete. 
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SUMMARY     AND 
RECOMMENDATIONS    FOR    FURTHER     RESEARCH 


SUMMARY 

This  thesis  constitutes  a  first  exploration  of  the 
use  of  interactive  graphics  for  conceptual  design  of 
submarines.  The  resulting  product  is  useful,  flexible, 
and  accurate  for  conceptual  design  activity.  It  is  by  no 
means  a  finished  product  in  the  commercial  sense,  but 
provides  a  tool  to  build  on  and  to  use  confidently  for 
such  applications  as  student  design  projects. 

The  power  of  real-time  graphic  display  of  the 
design  during  its  creation  is  a  major  factor  in  the 
worth  and  potential  of  the  program.  Many  hours  of 
laborious    manual    drafting    are  avoided    by    the 

high-resolution  graphic  output.  In  addition,  a  user  with 
even  a  rudimentary  knowledge  of  the  Computervision 
system  can  go  far  beyond  the  capabilities  of  the  basic 
software  package,  thus  greatly  enhancing  its  value  as  a 
design  tool . 
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The  program  is  highly  interactive;  the  user 
controls  all  important  design  decisions.  The  designer 
must  have  a  solid  knowledge  of  naval  architecture,  which 
is  appropriate  for  any  ship  design  program.  Subjective 
design  Judgements  must  remain  in  the  province  of  the 
user  for  a  credible  result,  and  that  philosophy  has  been 
rigorously  adhered  to. 

RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

There  a.re  nearly  unlimited  opportunities  for 
continued  work  on  this  project,  and  in  this  general 
area.  This  thesis  has  scratched  the  surface  in  an  area 
ripe  for  research  and  development. 

The  first  obvious  areas  encompass  those  portions  of 
the  design  process  which  were  excluded  from  the  current 
project.  These  include  structures,  propeller  design  and 
detailed  interior  arrangements. 

Submarine  structural  design  could  benefit  greatly 
from  treatment  in  the  CAD  environment.  Visualization  of 
complex  structural  relationships  is  intrinsically 
clearer  with  the  use  of  real-time  graphic  display.  In 
addition,  the  use  of  powerful  finite  element  packages, 
with  automatic  mesh  generation,  is  possible  on  many  CAD 
systems,  including  Computervision .  These  packages 
display  stress  levels  in  color,  and  will  magnify 
displacements  on  command  for  rapid  visual  analysis. 
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Arrangement  of  interior  details  is  feasible  an 
these  systems,  including  automatic  interference  checking 
and  layout  of  piping  and  electrical  distribution 
networks. 

Propeller  design  may  be  facilitated  by  graphic 
display  of  the  flow  regime  during  analysis.  Much 
analytical  work  remains  in  this  ar&a-,  and  research  into 
links  with  main-frame  systems  and  large  data  bases  is 
also  needed. 

In  addition,  the  translation  of  this  package  into 
the  Fortran— S  language,  with  attendant  programming  of 
linking  and  loading  files,  would  greatly  improve 
execution  speed  and  I/O  flexibility.  The  creation  of 
customized  graphics  commands  would  add  more  power  and 
enhance  the  ergonomics  of  the  package. 

Further  work  in  the  subject  areas  of  weight 
estimation,  hydrodynamics,  and  animated  display  of 
control  system  performance  is  required  to  derive  the 
maximum  benefit  from  submarine  CAD.  All  of  these 
subjects  have  been  separately  addressed  in  other 
software  packages,  but  would  be  valuable  additions  to  an 
integrated  design  program. 

With  increased  interest  and  design  volume  in  the 
submarine  design  field,  research  funding  is  both 
necessary  and  clearly  justified  by  the  potential 
improvements  to  the  speed  and  accuracy  of  the  design 
process.  The  overall  area  of  computer   aided   design   of 
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submarines  can  support  major  research  projects  in  the 
future,  benefitting  both  academic  and  commercial  design 
act i vi  t  ies. 
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APPENDIX  I : 
SAMPLE  OUTPUT 
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M IT. &BCD. WEIGHT 
3-28-84  13«37i44  FUTIL 


6.  18 


I 

2 
3 

4 
5 
6 
7 
8 
9 
JO 
II 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 


WEIGHT  ESTIMATE  FOR  13.461  DESIGN  PROJECT* 

TYPE*    ATTACK 

PR0PULSIX)N«    MODIFIED  WANKEL   ELECTRIC 

BATTERY   TYPE!    NICKEL   CADMIUM 

SPRINT    KW-Hi    3300 

ENDURANCE    KH=H»    8000 

BATTERY   WEIGHT*    165.631 

ENGINE   HORSEPOWER*    4500 

ENGINE    4    ASSOC.    EQUIP.    WEIGHT*    22.86 

MOTOR    &    ASS-OC.    EX1UIP.    WEIGHT*    32.25 

DISCRETE    ADDITION    TO  GROUP    2    WEIGHT*    80 

MAX.    OPERATING   DEPTH    IN  FEET*    700 

LEAD   FRACTION!    0. 1 
VARIABLE    LOAD  FRACTION*    0.089 
RESERVE    BUOYANCY    FRACTION!    0..I25 
FREE    FLOOD  FRACTION!    0.08 


CATEGORY 


1 

2 
3 
4 
5 
6 
7 
Al 


GROUP 

GROUP 

GROUP 

GROUP 

GROUP 

GROUP 

GROUP 

COND. 

LEAD 

COND.    A 

VAR.    LOAD 

NSC 

MBT 

SUB.    DISP. 

FREE    FLD. 

ENV.    DISP. 


WEIGHT 

493.318 
300.741 
17.4418 
95 

137.043 
I  17.3 
85 

1245.84 
124.584 
1370.43 
133.884 
J504.31 
1 88.039 
1692.35 
135.388 
1827.74 
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MlT.&BCD.GEOMOUT 
3-28-84    !3-«38«56   FUTIL 


6.  18 
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19 


GEOMETRY   OUTPUT   FOR    13.461    DESIGN    PROJECT! 


LENGTH    OVERALL!     139.947 

MID-BODY    LENGTH'     74.9465 

OVERALL   PRISMATIC!    0.798368 

LENGTH/DlAMETERi    5.1832 

SUBMERGED   LCB    (AFT   OF   FP )  «    66.2489 

DRAFT    0  .NORMAL   SURF.    COND.  »    21.9585 

LCB    9    NORMAL   SURF.    COND.       *    66.3329 

RECAP  OF  INPUT  DATA! 

ENVELOPE  DISPLACEMENT  1 1  827. 74 

ENTRANCE:  25 

RUN      i  40 

DIAMETER*  27 

FWD  PRISMATIC«  0.6 

AFT  PRISMATIC*  0.55 
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MIT.&BCD.OFFOUT 
3-28-84    13  08*56    FUTIL 


6.  18 


1! LENGTH 

2!0 

3 

.0.5 

4 

M.5 

5 

2.< 

59946 

6 

4.< 

)9839 

7 

6.99732 

8 

13 

9946 

9 

20 

992 

10 

27. 

9893 

1  1 

34. 

9866 

12 

41. 

984 

13 

48. 

9813 

141 

55, 

9  786 

15 

62 

,9759 

J6 

.69. 

9733 

17! 

76 

9  706 

18 

83 

9679 

.19 

90. 

9652 

20. 

97. 

9625 

21 

104.96 

22 

1  1 

.957 

23. 

118.955 

24 

125.952 

25 

132.949 

26! 

139.947 

HRT   FP 


OFFSET 


1.79298 
3.41525 

5.08919 

6.77129 

8. 1157 

II  .36/4 
13. 1284 

13.5 

13.5 
13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 
I3.J34 

12.444 

10.7078 

8.07667 

4.51632 

0 
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M IT. &BCD. POWER 
4-  9-84  10i13«59 


FUTIL   6.18 


1 

! SPEED   AND  J 

3X)WER   RESULTS    F 

OR    13.461 

2 

! 

3 

! NON-PROPULSION    LOADS    (KW) 

t     135 

4 
5 

! PROPULSIVE 

COEFFICIENT 

i    0.77 

6 
7 

8 
9 

SPEED    <KT> 

SHP 

TOTAL    KW 

1 

0.615861 

135.459 

10 

!    2 

4.61783 

138.444 

1  1 

!    3 

15.0496 

146.222 

12 

4 

34.8408 

160.981 

13 

5 

66.8572 

184.855 

14 

!    6 

113.92 

219.95 

15 

7 

178.818 

268.345 

16 

8 

264.312 

332.097 

17 

9 

373.14  1 

413.251 

18 

10 

508.023 

513.833 

19 

II 

67 1 . 663 

635.859 

20 

12 

866.75 

78J .335 

21 

13 

1095.96 

952.256 

22 

14 

1361.95 

1  150.61 

23 

15 

1667.38 

.1378.37 

24 

16 

2014.9 

1637.51 

25 

17 

2407.  13 

1930 

26 

18 

2846.71 

2257.79 

27 

19 

3336.24 

2622.83 

28 

20 

3878.34 

3027.08 

29 

21 

4475.63 

3472.47 

30 

22 

5130.68 

3960.95 

31  . 

23 

5846.  1  1 

4494. 44 

32 

24 

6624.47 

5074.87 

33! 

25 

7468.37 

5704.  17 

34 

26 

8380.38 

6384.25 

35 

27 

9363.07 

71  17.04 

36! 

28 

10419 

7904.45 

37! 

29 

1 1550.7 

8748.37 

38 

30 

12760.8 

9650.74 

39! 

31 

1405  1.8 

10613.4 

40 

32 

15426.3 

1  1638.4 

41 

33 

16886.8 

12727.5 

42 

34 

18435.8 

13882.6 

43 

35 

20075.9 

15105.6 

44 

36 

24  809.6 

16398.4 

45! 

37 

23639.5 

1 7762.9 

46' 

38 

25568 

19201 . 1 

47 

39 

2759  7.7 

207.14.6 

48! 

40 

29  731.2 

22305.5 
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MIT.&BCD.BALOUT 

5-   6-84    IO«54i3/    FUTIL   6.21 


1 

2 
3 
4 
5 
0 
/ 

a 

o 

10 

1 1 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
20 


RESULTS  FROM  BALANCE  MODULE 
GROUP   WEIGHT      LCG      VCG 


2 

3 
4 
5 
6 
7 
Al 


860.354 

623.431 

99.d|4 

105 

243.99 

195.5 

90 

22  lb. 09 


LEAD    221 .d09 


A 

V    LD 

MBT 


2439.9 
183. 64d 
327.943 


I  10 
120 
I  12 
60 
I  10 
/O 
90 

106. 19 
76 . dd9 
103.53 
90 
101.39 


15.5 

10 

15 

20.5 

9.5 

22.5 

10 

13.902 

8.1414 

13.378 

16 

15.5 


MARGIN    LEAD    (TONS)*    50    GVCG=D/2 .    LCG=L()A/2 
STABILITY    LEA)    (TONS):     171.809    (*VCG=    6    FT 
STABILITY    LEAD    REQUIRED    LCG*    67.1016 

NSC    (TONS>»    2623.55    ^LCB*     103.145 

SUBM  <TONS)«  2951.49  ;<*LCB«  102.951 

SUBMERGED   STABILITY    (BG):     1.72272    FT 
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MIT.OUT.&BCD.POLY 
5-  7_84  I9t24:02  FUTIL  6, 


I 

2 

3 

4 

5 

6 

/ 

8 

9 

10 

1  1 

12 

13 

14 

15 

16 

17 

lb 

19 

20 

21 

22 

23 

24 

25 

26 

27 

26 

29 

30 

31 


POLYGON  SUMMARY  OUTPUT  FILE 


ITEM 

CREW  AND  EFFECTS 
SLBM'S  OR  COMP.  WATER 
SANITARY  TANKS 
LUBE  OIL  IN  SUMPS 
FIXED  CLEAN  FUEL  OIL 
♦♦♦♦TOTAL  FIXED  ITEMS^*^* 
PROVISIONS  AND  STORES 
REVITALIZATION  02 
TORPEDOS  IN  FWD  RfX)M 
TORPEDOS  IN  AFT  R(X)M 
TACTICAL  MISSILES  FWD 
TACTICAL  MISSILES  AFT 
WRT  TANKS 

RESERVE  ELECTROLYTE 
TOTAL  FRESH  WATER 
RESERVE  LUBE  OIL 
FUEL  BALLAST  FWD 
FUEL  BALLAST  AFT 
VARIABLE  FUEL  OIL 


VARIABLE  BALLAST  TANK  DAT> 

TANK  GROUP 
FORWARD  TRIM 
AUXILIARY 
AFTER  TRIM 


LCG 


WEIGHT 


55 

6 

0 

0 

85 

1.5 

90 

5 

0 

0 

78.913 

1  1. 

75 

9.5 

88 

1 

45 

21 

0 

0 

0 

0 

1  16 

13 

30 

5 

58 

1 

100 

1  1 

95 

5.5 

26.5 

28 

1  15 

21 .074 

45 

49.074 

* 

LCG 

CAPA 

35 

25 

75 

10 

105 

35 
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APPENDIX  II: 
SOURCE  CODE  LISTINGS 


-16- 


MIT.^BCD. DESIGN 
5-    7-«4    13:50:00    FUTIL   0.21 


5HT    (TONS). 


FEET:     )ZMAX 


(3)  PROPULSION? 


DIM  W(l7),P(l7),F(d).RS(2l),ClRC(2l),LX<8).RX(d).XPOS<30).SA(2l)<# 

>,LS(2I  )  .HP (4  I  )  . KW(4  1  )  ,I)U(2I  )  ,  X(  29  )  .  Y (  29 ) 

READ  (ENTER  DESIGN  NAME:  )  *NAME 

^WEIGHT  CONTINUE 

N=l 

REPEAT 

N=N+I 

IF  (N.EQ.3.0R.N.EQ.5)  GOTO  I_fX)P 

PRNT  DO  YOU  WANT  TO  SPECIFY  (S)  OR  CALCULATE  (C>  GROUP  IN]  WEIGHTS? 

READ  *R 

IF  (&R.FQ."S")  F(N)=1 

IF  <*R.EQ."C")  F(N)=0 

IF  (F(N).EQ.O)  GOTO  LOOP 

PRNT  ENTER  GROUP  IN]  WEK 

READ  WT 

W(N)=WT 
#L(X)P  CONTINUE 
UNTIL(N.EQ.7) 

READCMX.  OPERATING  OEPTH  IN 
P( 1 )=.2I55+1.606  225E-4*ZMAX 
READ(RESERVE  BUOYANCY  FRACTION:  )P( 13) 
REAO(LEAD  FRACTION*  )P(Q) 
READ(VARIABLE  LOAD  FRACTION:  )P(I|) 
REA^<FREE  FLOOD  FRACTION:  )P(15) 
PRNT  NUCLEAR  (1).  DIESEL  (2).  OR  WANKEL 
READ  PM 

PRNT  ATTACK  (1)  OR  FBM  (2)7 
READ  T 

IF  (F(2) .EQ. I )  GOTO  LABI 
GOSUB  GR0UP2 

*LABI  IF  (F(6) .EQ. 1 )  GOTO  Al 
GOSUB  GR0UP6 
GOTO  Al 
#GR0UP2 

READ  (SHAFT  HORSEPOWER:  )SHP 

IF  (PM.NE. 1 )  GOTO  BATT 

K=.03*SHP/I000+I. 13 

WI=(K*1E5*SHP)/(LG(SHP ) ) **b 

W(2)=WI/2240 
RTNSUB 
#BATT  CONTINUE 

PRNT  BATrERY  TYPE:  NICAD  (1). 

READ  BT 

READ(T0TAL  ENGINE  HORSEPOWER: 

IF  (PM.EQ.2)  KE=.0096d 

IF  (PM.EQ.3)  KE=.0050d 

WE=KE*ESP 

WM=.00645*SHR 

READ(T()TAL  KW-H  STORAGE  FOR 

READ (TOTAL  KW-H  STORAGE  FOR 

IF  (BT.EO. I )  CS=35.b6lO 


VARTA  (2),  OR  TRIDENT  (3)7 


)ESP 


SPRINT:  )KSP 
ENDURANCE  RATE: 


)KEND 


IF  (BT.FQ.2) 
IF  (BT.EO. 3) 
B=KSP/CS 
IF  (BT.EO. 1 ) 
IF  (BT.FQ.2) 


CS=22.3b3 
CS=I2. 1^2b 


Ct=4d, 
CE=4^, 


3 
674 
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57 
5b 
50 
60 
61 
62 
63 
04 
65 
66 
6/ 
6d 
69 
70 
71 
12 
73 
IA 
lb 
lb 
77 
Id 
70 
80 

a  i 

a2 
a  3 
d4 
d5 
d6 
61 
dd 
ao 
oo 
oi 

02 

03 

04 

05 

06 

07 

Od 

00 

100 

101 

102 

103 

104 

105 

106 

107 

lOd 

100 

I  10 

I  I  I 

1  12 

I  13 

1  14 

I  lb 

!  \6 


IF    (HT.EQ.3)    CE=2t>.4012 

BI=KEND/CE 

IF  (Bl  .GT.R)  B  =  BI 

WT=WE+WM+B 

PRNT  ENTER  DISCRETE  ADJUSTMENT  TO  GROUP  2  (0  OR  #  OF  TONS) 

REAU  DA 

W(2)=WT+DA 
RTNSUB 
#GR0UP6 

HEAIXNUMBER    OF    PERSONNELi     )    NP 
rt  (6)=2.3*NP 
RTNSUB 
#AI 

IF    (PM.EQ. I )    P(3)=.045 
IF    ( PM.EQ. 2. OR. PM.EQ. 3)    P(3)=.0I4 
P(4)=.06 

IF    (T.EQ. 1  )    P(5)=.  1  I 
IF    (T.EQ. 2)    P(5)=.0d5 
IF    (T.EO. 1 .AND. PM.EQ. 1 )    P(7)=.03 
IF    (T.EQ. I .AND. PM.EQ. 2)    P(/)=.04 
IF    (T.EQ. I .AND. PM.EQ. 3)    P(7)=.05 
IF    (T.EQ. 2)    P(  /)=.  14 

FW=(P(  I )+P(  I )*P(0)  )/(  l-P(  I  I  ))+P(3)+P(5> 
IF  <F(4) .EQ.O. AND. F<7) .EQ.O)  GOTOCASEl 
IF  (F(4) .EQ. I .AND. F(7) .EQ.O)  G0T0CASE2 
IF  (F(4)  .EO.O.  AND.F(7)  .EQ.  1  )  G(TT0CASE3 
IF  (F(4) .EQ. I . AND.F17) .EQ. I )  GOTO  CASE4 
*CASE1 

A  I  =  (W(^)+W(6 ))/ ( 1 -FW-P(4 )-P(  / ) ) 

W(4)=P(4)*A1 

W(  /)=P(  /)*AI 
GOTO  TOTAL 
#CASE2 

AI=(W(2)+W(4)+W(6))/( l-FW-P( 7  )  ) 
rt( 7)=P( 7)*A1 
GOTO  TOTAL 
*CASE3 

Al =(W(2)+W(6J+W (  /))/(! -FW-P(4  )) 

rt(4)=P(4)*Al 
GOTO  TOTAL 
#CASE4 

A1=(W(2)+W(4)+W(6)+W( 7 ))/( 1-FW) 
#T0TAL 
rf(3)=P(3)*Al 
W(5)=P(5)*AI 
W(d)=A l 
W(0)=P(0)*AI 
W(  10)=A1  -t-W(O) 
W(12)=W( 10)/( 1  — P ( 1 1  )) 
W(  I  I  )=P(  I  I  )*W(  12) 
rt( I )=P( I )*W(12) 
W( 13)=P( I  3) *W ( 12) 
W( I4)=W( I2)+W( 13) 
ft  (  15)=P(  I5)*W(  14  ) 
W( I6)=W( 15)+W( 14) 
PRNT  CALCULATIONS  COMPLETE... 
PRNT 
PRNT 

<»»  OUTPUT  SECTION 
Ic  (T.EQ.  I  )  <T=MATTACK" 
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IF    (T.E0.2)    &T="BALLISTIC    MISSILE" 

IF    (PM.EO.I)    &PM=»NUCLEAR" 

IF  (PM.EQ.2)  &PM="DIESEL  ELECTRIC" 

IF  (PM.EQ.3)  &PM="WANKEL  ELECTRIC" 

IF  (DA.NE.O)  &M="M()OIFIED  " 

IF  (DA.EQ.O)  &M="" 

IF  (Bt'.EQ*.  I)  &BT="N1CKEL  CADMIUM" 

IF  (BT.EQ.2)  &BT="LEAD  ACID  (VARTA)" 

IF  (3T.EQ.3)  &BT="LEAD  ACID  (TRIDENT)" 

PHNT  WEIGHT  ESTIMATE  FOR  [&NAMEI* 

PRNT 

PRNT  TYPE*  [Ml 

PRNT  PROPULSION*  [&MH&PM] 

IF  (PM.EQ. I )  GOTO  SHAFT 

PRNT  BATTERY  TYPE*  UBTJ 

PRNT  SPRINT  KW-H*  (KSPJ 

PRNT  ENDURANCE  KW-H*  [KEND] 

PRNT  BATTERY  WEIGHT*  IB! 

PRNT  ENGINE  HORSEPOWER*  [ ESP  J 

PRNT  ENGINE  WEIGHT*  [WEI 

PRNT  MOTOR  WEIGHT*  [WMJ 

IF  (DA.NE.O)  PRNT  DISCRETE  ADDITION  TO  GROUP  2*  IDA] 

#SHAFT  CONTINUE 

PRNT  SHAFT  HORSEPOWER*  ISHP) 

PRNT  MAX.  OPERATING  DEPTH  <FT> *  CZMAX] 

IF  (F(6).EQ.I)  PRNT  NUMBER  OF  PERSONNEL*  [NPJ 

PRNT  LEAD  FRACTION*  rP(O)] 

PRNT  VARIABLE  LOAD  FRACTION*  (P(ll)] 

PRNT  RESERVE  BUOYANCY  FRACTION*  [P(I3)I 

PRNT  FREE  FLOOD  FRACTION*  [P(I5)I 

PRNT  HIT  <RETURN>  TO  CONTINUE... 

READ  ACONT 

PRNT 

PRNT  CATEGORY     WEIGHT 

PRNT 

1=1 

REPEAT 

PRNT  GROUP  [ I ]    (W(  I)  ] 

1=1  +  1 

UNTIL  (  I.EQ.d) 

PRNT  COND.  Al     [W(d) J 

PRNT  LEAD        IW(9) ] 

PRNT  COND.  A      [W(  10) ] 

PRNT  VAR.  LOAD    fW(  I  I  ) ] 

PRNT  NSC  [W(12) ] 

PRNT  MBT  [W( 13) I 

PRNT  SUB.  DISP.   [W( 14) ] 

PRNT  FREE  FLD.    [ W( 15) J 

PRNT  ENV.  DISP.   [W(I0)J 

PRNT 

PRNT  SATISFACTORY  (S)  Ok  RECALCULATE  (R>? 

READ  &DECID 

IF  (ADECID.EQ. "R")  GOTO  WEIGHT 

PRNT  PLEASE  WAIT  WHILE  OUTPUT  AND  PASS  FILES  ARF  WRITTEN.  THE  OUTPUT  FOR 

PRNT  THF  WEIGHT  CALCULATIONS  WILL  BE  IN  FILE  "M IT. WEIGHT" .  WHEN  YOU  HAVE 

PRNT  COMPLETE!!  THE  ENTIRE  DESIGN  PROGRAM.  "MIT. WEIGHT"  AND  THE  OUTPUT 

PRNT  FILES  FOR  SUBSEQUENT  MODULES  MAY  BE  OB  I'M  NED  IN  HARD  COPY  BY  USING 

PRNT  THE  COMMAND  'PRI NTLP"FI LFNAME"' . 

OPENW  2."MIT.WHIGHT" 


-ai 


I  77 
17b 
I7Q 
180 
Itil 
182 
183 
184 
l«5 
Ib6 
187 
I  88 
Iti9 
190 
191 
192 
193 
1Q4 
195 
196 
197 
196 
|9Q 
200 
201 
202 
203 
204 
205 
206 
20/ 
20o 
209 
210 

21  I 
212 
213 
214 
215 
216 
217 
218 
219 
220 

22  1 
??2 
223 
224 
225 
226 
227 
?2a 
220 
230 
23) 
232 
233 
234 
235 
P  36 


«.SHP=SHP 

&ESP=ESP 

4KSP=KSP 

&KEND=KEND 

&ZMAX=ZMAX 

&NP=NP 

&BB=B 

*WM=WM 

&WE=WE 

&DA=DA 

&LF=P(9) 

*VLF=P(  I  I  ) 

&RBF=P( 13) 

&FFF=P(I5) 

4WI=W( I ) 

&W2=W(2) 

&W3=W(3) 

&W4=W<4) 

&W5=W(5) 

*W6=W(6) 

&W7=W( 7) 

&W8=W (8) 

&W4=W(9) 

&wio=W( 10) 

&W I  | =W  <  I  I  ) 

4WI2=W( 12) 
&WI3=W( 13) 
&W1 4=W( 14) 
&WI5=W( 15) 
&W16=W( 16) 

&X="WEICHT  ESTIMATE  FOR  "♦&NAME+ 

GOSUB  0 

WRITEF  2,&P 

4X="TYPE:  "+&T 

GOSUB  0 

*X="PR0PULS10N:  "+&M+APM 

IF  (PM.EQ. 1 )  GOTO  NUKE 

GOSUB  0 

&X="BATTERY  TYPE:  "+4BT 

GOSUB  O 

&X="SPRINT  KW-H:  "+&KSP 

GOSUB  O 

AX="ENDURANCE  KW=H:  "+4KEND 

GOSUB  0 

*X  =  ''BATTERY  WEIGHT*  "  +  &BB 

GOSUB  0 

iX=»ENGINE  HORSEP()WER«  "+&ESP 

GOSUB  O 

&X="ENGINE    &    ASSOC.    EQUIP.    WEIGHT*    "+&WE 

GOSUB    0 

4X="M()T0R  4  ASSOC.  EQUIP.  WEIGHT'  "  +  &WM 

GOSUB  0 

IF  (OA.EQ.J)  GOTO  NUKE 

&X="DISCHETF.  ADDITION  TO  GROUP  2  WEIGHT'  "+*DA 

GOSUB  O 

#NUKE  CONTINUE 

4X=»MAX.  OPERATING  DEPTH  IN  FEET:  "+4ZMAX 

GOSUB  O 

WRITEF  2,*P 

RX=»I.FAD  FRACTION:  "+\LF 
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237JGOSUB  0 

23dJAX="VARIABLE  LOAD  FRACTION:  "+AVLF 

239 JGOSUB  0 

240!&X="RESERVE  BUOYANCY  FRACTION*  "+ARBF 

24  1 JGOSUB  0 

242!AX="FREE 

FL(X)D 

FRACTION*  "+AFFF 

243!G()SUB  0 

244IWRITEF  2, 

&P 

245JWRITEF  2, 

AP 

24o!&X="CATEGORY 

WEIGHT" 

24  /JGOSUB  0 

24d!WRITEF  2. 

AP 

249!AX="GROUF 

1 

"+AWI 

250 JGOSUB  0 

25! !AX=»GROUF 

2 

"+AW2 

252!GOSUB  0 

253!&X="GROUP 

3 

"+AW3 

254!G()SUB  0 

255!AX="GROUF 

4 

"+AW4 

256 JGOSUB  0 

257!AX="GROUP 

5 

"+AW5 

25d!GOSUB  0 

259!&X="GROUP 

6 

"+AW6 

260JGOSUB  0 

26  1 !AX="GROUP 

7 

"  +  &W7 

262!GOSUB  0 

263!AX="COND. 

Al 

"+AWd 

264 JGOSUB  0 

265!AX="LEAD 

"♦AWQ 

266 JGOSUB  0 

26/!&X="COND. 

A 

"+AWIO 

268 JGOSUB  0 

269!&X="VAH. 

LOAD 

"+AW I  I 

270JGOSUB  0 

271  !AX="NSC 

"+AWI2 

272 JGOSUB  0 

273!&X="M3T 

"+AWI3 

2/4 JGOSUB  0 

275!&X="5UB. 

OISP. 

"+AWI4 

2/6  JGOSUB  0 

277!AX="FREE 

FLD. 

"+AWI 5 

2 /d JGOSUB  0 

279!AX="ENV. 

DISP. 

"+AWI6 

2dO!GOSUB  0 

281 JGOTO  DONE 

2d2!#0  CONTINUE 

2d3JWRITEF  2, 

AX 

2d4!RTNSUB 

2d5!#D()NE  CONTINUE 

2ti6!OPENW  3," 

MIT.PASSI " 

2d/JWHITEF  3. 

AW  12 

2Hd!WRITEF  3. 

&WI5 

2dU!.1RITEF  3. 

AWIo 

2qojwri;ef  3. 

*NAME 

2^1  JOPENW  4," 

MIT. BALPASS" 

242JWRITEF  4, 

*W| 

2Q3JWRITEF  4, 

AW? 

2<J4!WRITEF  4, 

AW  3 

245JWRITEF  4. 

AW4 

296!WRITFF  4. 

K/ih 
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29  7 

WRITEF  4.  &W6 

29d 

WRITEF  4,  KHI 

299 

WRITEF  4,  &Wd 

300 

WRITEF  4,  &W9 

301 

.WRITEF  4.  &W10 

302 

WRITEF  4,  &W| I 

303 

WRITEF  4.  AWI2 

304 

WRITEF  4,  &WI3 

305 

WRITEF  4,  A.WI4 

306 

WRITEF  4,  &W15 

30/ 

WRITEF  4,  &WI6 

30d 

OPENW  5,"MIT.PASSPH2" 

309! 

WRITEF  5.&WI0 

310 

WRITEF  5.&NAME 

311. 

RUN  NEW  MIT.G3 

312 

END 
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MIT.&BCD.G3 
4-27-84  11*08*26  FUTIL  6.21 


I 
2 

3 

4 
5 
6 
7 
8 
Q 
10 
1  1 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 


OF  ENTRANCE*  )LF 
LENGTH  OF  HUN:  )LA 
HULL  DIAMETER*  )D 
FOREBODY  PRISMATIC  COEFF.* 
AFTERBODY  PRISMATIC  COEFF. 


DIM  RS(2I ), CIRC (2 1 ) ,LX (8) , RX (8 ) , XPOSC30 ) , SA( 21  )  ,LS(2I )<# 
>,HP(41 ),KW(4  1 ),DU(2I ) ,Xr I32).Y( 132) ,SX(21 ) 

#GEOM 

OPENR  I  , "MIT. PASS  I" 

READF  1,&W 

DNSC=&W 

READF  l.&W 

DFF=&W 

READF  1,&W 

DENV=&W 

READF  l.&NAME 

PI=3.  14159265 

#BRANCH  CONTINUE 

READ(ENTER  LENGTH 

READ( ENTER 

READ( ENTER 

READ( ENTER  FOREBODY  PRISMATIC  COEFF. *  )CPF 

READ( ENTER  AFTERBODY  PRISMATIC  COEFF. »  )CPA 

L=LF+LA 

LL=LF 

lNTV=L/20 

R  =  D/2 

NF=71.04  77*CPF**5+65.8J07*CPF**4-483.3037*CPF**3+587.4I37*CPF**2<# 

>-28) .7224*CPF+4Q.5876 

NA=379.6546*CPA**5-938.4  708*CPA**4+9  44.8853*CPA**3-4  71 .0872*CPA**2<# 

>+l 19. |465*CPA-1 1 .3454 

LS( I )=0 

LS<2! )=L 

J=2 

REPEAT 

LS(J)=(J-1 )*INTV 

XL=LS(J) 

GOSUB  RAD 

RS< J)=RANS 

SA(  J)=PI*RS( J)**2 

DU(J)=SA(J) 

J=J+I 

UNTILU.EQ.2l  ) 

GOSUB    INTG 

VENV=35*DENV 

DV=VENV-ANS 

LMB=DV/(PI*R**2) 

LL=LF+LM3 

LOA=LF+LMB+LA 

INTV=LOA/20 

LS(21  )=L()A 

LS< I )=0 

J  =  2 

REPEAT 

LS(J)=(J-1 )*INTV 
XL=LS( J) 
GOSUB    RAD 

RS( J)=RANS 

CIRC( J)=2*PI*HS(J) 

SA(J)=PI*RS(J)**2 

DU(.D=CIRC(J) 
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57 

3d 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

g4 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

I  10 

I  I  I 

I  12 

I  13 

1  14 

1  15 

I  16 


J=J*I 

UNTIU  J.EQ.21  ) 

GOSUB    INTG 

SF=ANS 

VCYL=PI*L0A*R**2 

CP=VENV/VCYL 

J  =  2 

REPEAT 

DU(J)=LS(J)*SA(J) 

J=J  +  I 

UNTILf J.EQ.21 ) 

GOSUB  INTG 

LCBSUB=ANS/VENV 

DR=.5*R 

0N=DNSC*.8*0FF 

L0=DN-.005*DN 

HI=DN+.005*DN 

PR  NT  ONi  [DN] 

PRNT  L0«  [LO] 

PRNT  HI»  [HI] 

#LCBLP  CONTINUE 

PRNT  ITERATING  ON  NSC  DRAFT  »  T=  [R+DRJ 

J=2 

REPEAT 

IF(DR.GE.RS(J  ))  GOTO  SAME 

AC=(AC()S(DR/RS(J  )))*(  PI/180) 

SX  (J)=SA(J)-(RS(J)**2*AC-DR* (RS ( J ) **2-DR**2) **. 5 ) 

GOTO  SKIP 

#SAME  CONTINUE 

SX( J)=SA(J) 

#SKIP    CONTINUE 

DU(J)=SX(J) 

J=J+1 

UNTIL(J.EQ.2I ) 

GOSUB  INTG 

DM=ANS/35 

IF(0M.GE.LO.AN0.OM.LE.HI )  GOTO  MOM 

IF(DM.GT.HI)  GOTO  MINUS 

DR=DR+. 1*DR 

GOTO  LCBLP 

#MINUS  CONTINUE 

DR=OR-. I*DR 

GOTO  LCBLP 

#MOM  CONTINUE 

J  =  2 

REPEAT 

DU(J)=LS(J)*SX (J) 

J  =  J  +  I 

UNTILU.EQ.21  ) 

GOSUB    INTG 

LCBNSC=ANS/(0M*35) 

UX=(INTV-2)/5 

LX( I )=.5 

LX (2) =1 

LX(3)=1 .5 

LX(4)=2 

REPEAT    LfX)P2s  J  =  5,  1  .  (J.GT.b) 

LX( J)=(J-4)*0X+2 

#LfX)P2   CONTINUE 

WEPEAT    EXTRA:  1=1  ,1  .(J.GT.b) 
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XL=LX(J> 

GOSUB    RAD 

RX( J)=RANS 

#EXTRA 

GOTO  GOUT 

#RAD    CONTINUE 

IF(XL.LE.LF)    GOTO   FRONT 

IF( XL.GT.LF.AND.XL.LE.LL)    RANS=R 

IF(XL.GT.LL)    GOTO   AFT 

GOTO    L(X)P4 

#FRoNT   CONTINUE 

XX=LF-XL 

RANS=R*(( 1-(  XX/LF)**NF)  **( l/NF)) 

GOTO    L(X)P4 

#AFT    CONTINUE 

XX=XL-LL 

RANS=R*( l-( XX/LA) **NA) 

#L(X)P4    RTNSUB 

#INTG    <#SIMPS0N'S    RULE    INTEGRATION 

ODD=0 

REPEAT   ODDL«J=2.2. (J.GT.20) 

ODD=onD+4*DU( J) 

#ODOL    CONTINUE 

EVEN=0 

REPEAT  EVENLsJ=3,2. (J.GT. IQ) 

EVEN=EVEN+2*DU(J ) 

*EVENL  CONTINUE 

ANS=( INTV/3)*(DU(  I )  «-ODD+EVEN+DU(  21  )) 

RTNSUB 

#GOUT 

LD=LOA/D 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

READ 


[DENV] 


ENVELOPE  DISPLACEMENT' 

LENGTH  OVERALL:  t LOA  ] 

OVERALL  PRISMATIC:  ICPJ 

DIAMETER*  [DI 

LENGTH/DIAMETER:  [LDI 

MIDBODY  LENGTH:  [LMB] 

FWD  PRISMATIC:  CCPF] 

AFT  PRISMATIC: 

ENTRANCE:  (LFJ 

RUN:  [LA] 

LCB  SUBMERGED: 

LCB  O  NSC     : 

DRAFT  9  NSC   : 

SATISFACTORY  (S)  OR  RECALCULATE  (R)? 

*DECID 

IF  (&DECID.EQ. "R")  GOTO  BRANCH 

PRNT    PLEASE    WAIT  WHILE    OUTPUT    AND    PASS    FILES    ARE    WRITTEN.    THE    OUTPUT 
PRNT    FILE    FOR    GEOMETRY    SUMMARY    DATA    WILL   BE    "M IT. GEOMOUT".     A    FILE    OF 
PRNT   OFFSETS    WILL   BE    IN    "M IT.OFFOUT" . 
OPENrt    3, »M IT.OFFOUT" 

&Q="LENCTH  WRT  FP  OFFSET" 

GOSUB  W 
LS(  I  )=U 
RS( I )=0 
&X=LS( I ) 
&C=RS< I ) 


fCPA) 


[LCBSUB] 

:  [LCBNSC] 

:  [R+DR] 

(S)  OR  RECALCULATE 
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177 
178 
179 

lao 

181 
182 
183 
184 
185 
186 
187 
188 
I  89 
190 
191 
192 
193 
194 
195 
196 
19/ 
198 
199 
200 
201 
202 
203 
204 
205 
206 

20  7 
208 
209 
210 

21  I 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 

22  7 
228 
229 
230 
231 
232 
233 
234 
235 
236 


&Q=AX+AS+AC 

WRITEF  3,&Q 

REPEAT  L(X)P5:J=I  .1  .(J.GT.8) 

&X=LX( J) 

&C=RX(J) 

4Q=&X+&S+&C 

GOSUB    W 

J=J+1 

#L(X)P5   CONTINUE 

REPEAT    LfX)P6«J=2.1  .(J.GT.2I  ) 

&X=LS(J) 

&C=RS( J) 

&Q=&X+AS+&C 

GOSUB    W 

#L(X)P6    CONTINUE 

GOTO    FILE 

#W    CONTINUE 

WRITEF   3.&Q 

RTNSUB 

#FILE    CONTINUE 

OPENW    4, "MIT. PASSHP" 

AX=CP 

WRITEF  4.&X 

AX=L()A 

WRITEF  4,&X 

AX=D 

WRITEF  4,&X 

AX  =  SF 

WRITEF  4.&X 

WRITEF  4.&NAME 

OPENW  5.J,MIT.GE()M0UT" 

AX="GEC)METRY  OUTPUT  FOR  "+&NAME+"»" 

WRITEF  5. AX 

WRITEF  5.&S 

WRITEF  5. AS 

&P=LOA 

&X="LENGTH  OVERALL*  "  +  &P 

WRITEF  5,&X 

&P=LMB 

&X="M ID-BODY  LENGTH:  "+AP 

WRITEF  5, AX 

AP=CP 

&X=-»OVERALL  PRISMATIC:  "+&P 

WRITEF  5. AX 

&P  =  LD 

&X="LENGTH/DIAMETER«  "+&P 

WRITEF  5, AX 

AP=LCBSUB 

AX="SUBMERGED  LCB  (AFT  OF  FP ) «  "+AP 

WRITEF  5. AX 

TT=R+DR 

AP=TT 

AX="DRAFT  0  NORMAL  SURF.  COND. «  "+AP 

WRITEF  5,&X 

&P=LCBNSC 

AX  =  "[_CB  9  NORMAL  SURF.  COND.   «  "+&P 

WRITEF  5. AX 

WRITEF  5. AS 

AX="RECAP  OF  INPUT  DATA*" 

WRITEF  5. AX 


23  7! 

WRITEF   5,&S 

238! 

&P=DENV 

239! 

&X="ENVELoPE 

DISPLACEMENT    :"+&P 

240! 

WRITEF    5,&X 

241! 

&P=LF 

242! 

&X=*'ENTRANCEt 

"+S.P 

24  3! 

WRITEF    5.&X 

244! 

&P=LA 

245! 

&X-"RUN               : 

"+&P 

246! 

WRITEF    5,&X 

247! 

&P=D 

24d! 

&X=J,nlAMETER: 

"+&P 

249! 

WRITEF    5.&X 

250! 

&P=CPF 

25  1  ! 

&X="FWD    PRISMATIC:    "+&P 

252! 

WRITEF    5.&X 

253! 

&P=CPA 

254! 

&X="AFT   PRISMATIC:     "+&P 

255! 

WRITEF    5.&X 

256! 

OPENW    6,J,MIT. 

LCBJ' 

257! 

&X=LCBNSC 

25d! 

WRITEF   6,&X 

259. 

&X=LCBSUB 

263! 

WRITEF   6. AX 

261 

&X=R 

262! 

WRITEF   6,&X 

263 

&X=LOA 

264! 

WRITEF    6.&X 

265 

OPENW    7. "MIT. 

PASSPHI" 

266! 

&x=n 

267 

WRITEF    7,&X 

26ti 

&X=LOA 

269 

WRITEF    7,&X 

2  70 

&X=LF 

271 

WRITEF    /,&X 

272 

*X  =  LA 

273. 

WRITEF    7.&X 

274 

&X=LCBSUB 

275 

WRITEF    7.&X 

2/6 

DSUB=DENV-DFF 

277 

&XOSUB 

2  76 

WRITEF    7.&X 

2  79 

.&X=DENV 

2d0 

WRITEF    /.&X 

281 

&X=NF 

262 

WRITEF    7,&X 

263 

!&X=NA 

284 

!WRITEF    7,&X 

265 

!#DRAW 

286 

!  X  (  1  ) =0 

28  7 

! Y ( 1 ) =0 

288 

!I=2 

289 

SREPEAT 

2«0 

!X( I)=LX( 1-1 ) 

29  1 

!Y( I)=RX( I-l ) 

292 

! 1=1+1 

293 

!UNTIL(  I.GT.9) 

294 

!LFINT=(LF-(INTV+I ))/23 

2Q5 

!XL=INTV+I 

?96 

!  1  =  10 
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297 

REPEAT 

29d 

X(I)=XL 

299 

GOSUB  RAD 

300 

Y(D=RANS 

301 

.XL=XL+LFINT 

302 

1=1  +  1 

303 

UNTIL( I.GT.3I ) 

304 

LINT=(L()A-LF)/101 

305 

XL=LF+LINT 

306 

1=32 

307 

REPEAT 

30d 

X( I)=XL 

309 

GOSUB  RAD 

310 

Y(I)=RANS 

31  1 

XL=XL+LINT 

312 

1=1  +  1 

313 

.UNTIL(  I.GT..I3I  ) 

314 

X(l  32)  =L()A 

315 

.Y( 132)=0 

316 

OPENW  d, "MIT. POINTS" 

317 

.J=l 

31d 

REPEAT 

319 

&X=X(J) 

320 

&Y=Y(J) 

321 

WRITEF  8,&X 

3  22 

WRITEF  8.&Y 

323 

J=J+1 

324 

UNTIL( J.GT. 132) 

325 

RUN  NEW  MIT. SPEED 

326 

END 
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MIT.&BCD. 
4-1  1-84 


SPEED 
I3i20«34 


FUTIL   6.18 


DIM  HP(40),KH(40) 

READt ENTER  HOTEL  AND  COMBAT  SYSTEM  LOADS  IN  KW» 

READ(ENTER  PROPULSIVE  COEFFICIENT!  ) PC 


)HL 


,-MIT.PASSHP" 

,&CP 

,&LOA 

.40 

,4SF 

,&NAME 


OPENR 

READF 

READF 

READF 

READF 

READF 

CP=&CP 

LOA=&L()A 

D=&D 

SF=&SF 

DA=I .09065E-3*LOA*D+I I .25 

J  =  l 

REPEAT 

RE  =  (  J*l  .68<?*L0A)/1  .27908E-5 

CF=.0  75/(LGfRE)-2)**2 

CR*CF*( 1 . 5**D/L0A)**I . 5+7*( D/LOA) **3) +. 002*(CP-.6) 

CT=CF+CR+. 00025 

HP(J)=( l/PC)*.00872*J**3*(CT*SF+DA> 

KH(J)=.7457*HP(J)+HL 

J=J+I 

UNTIL( J.GT.40) 

PRNT 

PRNT 

PRNT 

PRNT  SPEED  AND  POWER  OUTPUT  FOR  [&NAMEJ 

PRNT 

PRNT  NON  PROPULSION  LOADS  (KW)»  [HL1 

PRNT  PROPULSIVE  COEFFICIENT    «  [PC] 

PRNT 

PRNT  HIT    <RETURN>  TO      SEE    OUTPUT   VALUES... 

READ   &CONT 

PRNT 

PRNT   SPEED    (KT)  SHP  TOTAL   KW 

PRNT 

J=l 

REPEAT 

PRNT    [J]  CHP(J)]  [KW(J) ] 

J=J  +  I 

UNTIL(J.GT.20) 

PRNT   HIT    <RETURN>    TO   CONTINUE... 

READ    &CONT 

PRNT   SPEED    (KT) 


PRNT 

REPEAT 

PRNT    IJ] 

J=J+I 

UNTILU.GT.40) 

PRNT  PLEASE  WAIT 

PRNT  MODULE  WILL 


SHP 


CHP(J>] 


TOTAL  KW 


IKW(J) ] 


WHILE  THE  OUTPUT  FILE 
BE  IN  FILE  "MIT. POWER". 

OPENW  2, "MIT. POWER" 

&X="SPEED  AND  POWER  RESULTS  FOR  "♦&NAME+":" 

WRITEF  2.&X 

&PADI="         " 


IS  CREATED.  OUTPUT  FOR  THIS 


-91- 


57!&PAD2="      " 

581WRITEF  2.&PADI 

59!&HL=HL 

60!4PC»PC 

6.1  J&X  —  NON-PROPULSION 

LOADS  ( 

62!HRITEF  2,&X 

631&X— PROPULSIVE  COEFFICIENT 

645HRJTEF  2.&X 

65!HRITEF  2.&PADI 

66!&X="SPEED  <KT> 

SHP 

67!WRITEF  2.&PAD1 

68!HRITEF  2.&X 

6QIWRITEF  2.&PADI 

70!  J=.l 

71 'REPEAT 

72!&J-J 

73!&X»HP(J) 

74!&HP«&X( 1 ,8) 

75!&Y«KW(J) 

76!4KW=&Y( 1 ,8) 

77!&OUT«"  -+&J+&PADI+"    "+&H 

78!HRITEF  2  ,&OUT 

7Q!  J=J+.1 

80!UNTIL(J.GT.40) 

81 !RUN  NEW  MIT.ENV 

82'END 

(1CW)«    -+&HL 
t    •"♦&PC 


TOTAL   KW" 
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MIT.&BCD. BALANCE 
4-27-84  I0«47«27  FUTIL  6.21 


1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
I  I 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
26 
26 
21 
2d 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 


DIM  W( I6).LC(  16)  ,VC<  16)  ,LM( 16)  ,VM( 16) ,VY(  16) 

OPENR  I  .-"MIT.BALPASS" 

J*  I 

REPEAT 

READF  l.&X 

W(J)=&X 

J=J+I 

UNTILU.GT.  16) 

OPENR  2,',MIT.LC3J, 

READF  2.&X 

LCBNSC=&X 

READF  2.&X 

LCBSUB=&X 

READF  2,&X 

R=&X 

READF  2.&X 

L()A=&X 

PRNT  THE  INPUT  DATA  FOR  THIS  MODULE  AND  THE  POLYGON  MODULE  WILL  BE 

PROMPTED  IN  A  NEW  FORMAT.  WHEN  THE  PROGRAM  PROMPTS  FOR  AN  INPUT 
VALUE,  THE  CURRENT  VALUE  WILL  BE  DISPLAYED  TO  THE  RIGHT  OF  THE 
INPUT  PROMPT  MESSAGE.  THIS  WILL  APPEAR  AS  « 


THIS  WILL  FACILITATE  CHANGES  TO  INPUT  DATA 
ATION  THROUGH  THE  MODULE(S)  IS  REQUIRED. 
CONTINUE 


PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

#START 

LT=0 

VT=0 

J=l 

REPEAT 

PRNT 

PRNT  GROUP  [J]  LCG« 

READ(VALUE?)LC(J) 

PRNT  GROUP  (J  I  VCG: 

READ* VALUE?)VC(J) 

LM(J)=W(J)*LC(J) 

VM(J)=W(J)*VC(J) 

LT=LT+LM(J) 

VT=VT+VM(J) 

J=J+I 

UNTILU.GT.  7) 

LC(8)=LT/W(8) 

VC(8)=VT/W(8) 

LM(8)=LT 

VM(8)=VT 

PRNT 

READ(LCG  OF  VARIABLE 

READCVCG  OF  VARIABLE 

PRNT 

LM(  I  I  )=W( I  1  )*LC(  I  I ) 

VM(  I  I ) =W( 1  I  )*VC( I  1 ) 

PRNT  TOTAL  LEAD  IS  [W(9>] 

PRNT 

READ(AMOUNT  OF  LEAD  FOR  MARGIN 

SL=W(9)-ML 


INPUT    PROMPT*         =  CURRENT    VALUE' 


IF    MORE    THAN    ONE    ITER- 


L()AD?)LC(  I  I  ) 
L()AD?)VC(  I  1  ) 


TONS, 


IN   TONS?)ML 


57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

12 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

9  / 

98 

9<J 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

1  10 

I  I  I 

1  12 

1  13 

I  14 

1  15 

I  16 
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LSL=(W( 12)*LCBNSC-LM(8)-LM(  I  I  )-ML*( L0A/2 ) )/SL 

LC(9)=(SL*LSL+ML*(L0A/2))/W<9) 

PRNT  REQUIRED  STABILITY  LEAD  LCG  IS  [LSLJ. 

READCIS  THIS  FEASIBLE  <Y/N>?  )&R 

IFC&R.EQ."N")  GOTO  START 

LMC9)=W(9)*LC(9) 

VM(9)=SL*6+ML*R 

VCC9)«VM(9)/W(9) 

LM(10)=LM(8)+LM(9) 

VM( 10)=VM(8)+VM(9) 

LC( 10)=LM( 10)/W( 10) 

VC( I0)=VM( 10)/H( 10) 

LM( 12)=LM( IO)+LM(  1  1  ) 

VM(12)=VM( 10J+VMC 1  I  ) 

VC( 13) =R 

VM( I3)=R*WC 13) 

LCC 13)=CW< I4)*LCBSUB-LM( I2))/W( 13) 

PRNT 

PRNT  REOUIRED  MBT  LCG  IS  [LCC 13)]. 

READCIS  THIS  FEASIBLE  <Y/N>?  )&R 

IF(&R.EQ."N")  GOTO  START 

LMC 13)=WC J3)*LCC 13) 

LMCI4)=LM( 12)+LM( 13) 

LC( |4)=LM( 14)/W( 14) 

VM(  |4)=VM( I2)+VM(  13) 

VC( 14)=VM( I4)/W<  14) 

BG=R-VC( 14) 

PRNT 

PRNT  SUBMERGED  STABILITY  (BG)  IS  IBG]  FT. 

READCIS  THIS  ACCEPTABLE  <Y/N>?  )&R 

IFC&R.EG. "N")    GOTO  START 

PRNT  THE  OUTPUT  FOR  THIS  MODULE  WILL  BE  IN  FILE  "MIT. BALOUT" . 

J=l 

REPEAT 

VYCJ)=VC(J)-R 

J=J+1 

UNTILCJ.GT. 14) 

OPENH  3,  J,MIT.  BALOUT" 

&X="RESULTS  FROM  BALANCE  MODULE" 

WRITEF   3,&X 

&P I ="  " 

&P2="         " 

&P3="      " 

WRITEF   3.&PI 

&X="GROUP      WEIGHT  LCG  VCG" 

WRITEF    3,&X 

WRITEF    3,&P1 

J=l 

REPEAT 

&J=J 

&X=W( J) 

&W=&X( I .8) 

&X=LC(J) 

&LC=&X( I  ,6) 

iX=VC(J) 

&VC=&X( I  ,6) 

&OUT="  "+&J+"     "+&W+&P2+&LC+&P3+&VC 

WRITEF    3,&0UT 

J=J+I 

UNTILCJ.GT.  /) 
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I  17 

118 

1  19 

120 

12J 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

1  44 

145 

146 

14/ 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

I  66 

16/ 

168 

169 

1.70 

171 

I  12 

173 

I  74 

I  75 

I  It 


J  =  8 

&G="  A  I 
GOSUB  10 

J=9 

&G="  LEAD  Jl 

GOSUB  10 

J=IO 

&G="  A     " 

GOSUB  10 

J  =  1I 

&G="  V  LD  " 

GOSUB  10 

J=  I  3 

&G="  MBT  " 

GOSUB  10 

&NSC=W( 12) 

&SUB=W(I4) 

&LCBN=LCBNSC 

&LCBS=LCBSUB 

&BG=BG 

WRITEF   3.&PJ 

&SL=SL 

&ML=ML 

&()UT="MARGIN    LEAD    (TONS* «    "+&ML+"    «*VCG=D/2,    LCG=L()A/2" 

WRITEF   3.&0UT 

&OUT="STABILITY  LEAD  <TONS>»  "+&SL+"  0VCG=  6  FT" 

WRITEF   3.&0UT 

&LSL=LSL 

AOUT="STABILITY  LEAD  REQUIRED  LCG«  "+&LSL 

WRITEF  3.&0UT 

WRITEF  3.&PI 

&OUT=-"NSC  (TONS)  i  "♦&NSC+J'  8LCB«  "+&LCBN 

WRITEF  3.&0UT 

WRITEF  3.4PI 

&0UT="SU3M  (TONS) i  "+&SUB+"  0LCB«  "+&LCBS 

WRITEF   3.&0UT 

WRITEF   3,&PI 

&OUT*" SUBMERGED  STABILITY  (BG)«  "+&BC+"  FT" 

WRITEF  3.&0UT 

OPENW  4."MIT.PASSP0LY" 

&X=LCBSUB 

WRITEF  4,&X 

&X=LC( 1 1 ) 

WRITEF  4,&X 

&X=LC( 13) 

WRITEF  4,&X 

GOTO  PLOT 

#10  CONTINUE 

&X=W( J) 

&W=&X( I ,8) 

&X=LC(J) 

&LC=&X( I .6) 

&X=VC(J) 

&VC=*X( I  ,6) 

&OUT-&G+&W+4P2+&LC+&P3+&VC 

WRITEF    3.MHJT 

RTNSUB 

#PL()T 

SELECT  MODE  MODEL 

SELECT  LAYER  3 


)  77 
178 
179 
I  BO 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
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j=i 

REPEAT 

&J*J 

INSERT  P()INT«XfLC(J)  J  Y  [  VY(  J  )  ]  <CR> 

XP=LC(J)-H 

INSERT  TEXT  I&JJ  HGT  I .5« X[ XP ) Y[ VY ( J ) ] <CR> 

J=J+I 

UNTIL<J.GT.7) 

INSERT  P()INT»X[LC(9)]YIVY(9)  ]<CR> 

XP=LC(9)+1 

INSERT  TEXT  'LEAHK  HGT  1 .5 « XI XP ] Y[ VY(9 ) J <CR> 

INSERT  P()INT«X[LC(  13)  IYIVY(  13)  )<CR> 

XP=LC( 13)+l 

INSERT  TEXT  'MBT'  HGT  1 .5 i XC XP3 Y[ VY ( 1 3 ) ] <CR> 

RUN  NEW  MIT. POLY 

END 
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MIT.4BCD.P0LY 

5-  7-d4  |9:|4:02  FUTIL  6.21 


2 
3 
4 
5 
6 
/ 
d 
g 

10 

1 1 

12 

13 

I  4 

15 

16 

1  / 

Id 

19 

20 

21 

22 

23 

24 

25 

26 

27 

2d 

29 

30 

31 

32 

33 

34 

35 

36 

3/ 

3d 

39 

40 

41 

42 

4  3 

44 

45 

46 

47 

48 

4<J 

50 

51 

5? 

53 

54 

55 

56 


DiM  W< 16) ,XV( 10) ,YV< 10) ,XP(I0) ,YP( 10) ,XT( 

OPENR  I ."MIT.BALPASS" 

J=l 

REPEAT 

READF  I, XX 

W(  J)=AX 

J=J+1 

UNTILU.GT. 16) 

VLD=W(  I  I ) 

MBT=W( 13) 

DSUB=W( 14) 

OPENR  2,"MIT.PASSP()LY" 

READF  2. AX 

LSUB=&X 

HEADF  2.&X 

LVLD=&X 

READF  2.&X 

LMBT=&X 

#START  CONTINUE 

BAL=VLD 

PRNT   TOTAL   VARIABLE    LOAD    IS    [VLDl    TONS. 

PrtNT 

READ(WEIGHT    OF   CREW    AND    EFFECTS?)CE 

READ(LCG?)LCE 

BAL=BAL-CE 

GOSUB    STAT 

READMEIGHT    OF 

READ(LCG?)LMC 

BAL=BAL-MC 

GOSUB   STAT 

READ(WEIGHT   OF 

READ(LCG?)LST 

BAL=3AL-ST 

GOSUB    STAT 

READ(WRIGHT    OF 

READCLCG?) LLO 

BAL=BAL-L() 

GOSUB    STAT 

READfWEIGHT  OF 

READ  (  LOG?  )LCF() 

BAL  =  BAL-CF() 

GOSUB  STAT 

KEADfWEIGHT  OF 

READ(LCG?)LPS 

BAL=BAL-PS 

GOSUB  STAT 

READ (WEIGHT  OF 

REA0(LCC?)L02 

BAL=BAL-02 

GOSUB  STAT 

RFAD(WE1GHT  OF 

HEAU(LCG?)LTF 

BAL=8AL-TF 

GOSUB  STAT 

4EADIWEIGHT  OF 

3EAI)(I.CG?)LTA 


O) ,YT( 10) ,M(d) ,Y(d) 


sLBM'S  OR  COMPENSATING  WATER?)MC 


SANITARY  TANKS  AND  WATErt?)ST 


LUBE  OIL  IN  SUMPS?)LO 


FIXED  CLFAN  F.  ().  <IE.  SHIELD  TANK>?)CFO 


PROVISIONS  AND  STORES? ) PS 


REV ITALIZATION    0XYGEN?)02 


TORPEDOS    IN    FORWARD    HOOM?)TF 


TfJRPEOOS    IN    AFT    RfX)M?)TA 
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57 
58 
59 
00 
61 
62 
63 
64 
o5 
66 
67 
6a 
0Q 

ra 

7! 
12 
73 
74 
75 
10 
77 
7b 

/q 
ao 
til 
d2 
d3 
d4 
85 
86 
87 
dd 
8Q 
90 
Q1 
92 
93 
y4 
^5 

4/ 
48 
Q9 
100 
101 
102 
103 
104 
105 
106 
10/ 
108 
104 
I  10 
I  I  I 
1  12 
I  13 
I  14 
1  15 
I  16 


BAL=BAL-TA 

GOSUB   STAT 

READ(WEIGHT  OF  TACTICAL  MISSILES  F()RWAHD?)MF 

READ(LCG?)LMF 

BAL=BAL-MF 

GOSUB  STAT 

READ(WEIGHT  OF  TACTICAL  MISSILES  AFT?>MA 

READ(LCG?)LMA 

BAL=3AL-MA 

GOSUB  STAT 

TM=TF+TA+MF+MA 

LTM  =  (TF*LTF+.TA*LTA+MF*LMF+MA*LMA)/TM 

READ( WEIGHT   OF   WATER    FOR    WRT    TANKS? )WRT 

READ(LCG?)LWRT 

BAL=BAL-WRT 

GOSUB    STAT 

READ(WEIGHT    OF    RESERVE    ELECTROLYTE?) EL 

READ<LCG?)LEL 

BAL=BAL-EL 

GOSUB    STAT 

READ(WEIGHT   OF    FRESH    WATER    <P()TTABLE  .FEED,  BATTERY    RESERVE>?)FW 

READ(LCG?)LFW 

BAL-BAL-FW 

GOSUB    STAT 

READ(WEIGHT  OF  RESERVE  LUBE  &  HYDRAULIC  OIL?)RL() 

READ(LCO?)LRL() 

BAL=BAL-RLO 

GOSUB  STAT 

PHNT  THE  LAST  VARIABLE  LOAD  ITEM  TO  BE  CONSIDERED  IS  FUEL  OIL.  TO  ENSURE 
A  MARGIN  OF  SAFETY  AT  THE  BASE  OF  THE  POLYGON,  A  VARIABLE  FUEL  OIL 
(VFO)  TANK  IS  AUTOMATICALLY  SPECIFIED.  THE  REQUIRED  VFO  CAPACITY  IS 
COMPUTED  USING  YOUR  INPUT  OF  THE  FRACTION  OF  TOTAL  FUEL  TO 
PLACE  INTO  THE  VFO  TANK.  THE  MINIMUM  SUGGESTED  FRACTION  IS 
.23  (THE  RATIO  DIFFERENTIAL  OF  FUEL  AND  WATER  DENSITIES). 
IT  IS  HIGHLY  RECOMMENDED  THAT  AN  ADDITIONAL  MARGIN  OF  AT 
LEAST  .05  BE  SPECIFIED  TO  ALLOW  A  SAFETY  MARGIN  AT  THE  BASE 


OF  FUEL  TO  6E    PLACED  IN  VFO?)  VFR 


PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT  OF  THE  POLYGON 

PRNT 

READ( FRACTION 

TFO=BAL 

VFO=VFR*TF() 

FBT=TFO-VF() 

FL=CE+MC+ST/3+L0+CF0 

LFL=( CE*LCE+MC*LMC+ST*LST+LO*LLO+CFo*LCFO) /FL 

PRNT  FBT  CAPACITY  (TONS)*  IFBTI 

PRNT  VFO  CAPACITY  (TONS): 

PRNT 

READ(  AMOUNT  OF  F3T  <T()NS> 

FBTA=FBT-FBTF 

PRNT 

READ(FWD  FBT  LCG?)LFHTF 

READ (AFT  F3T  LCG?)LF3TA 

LFKT=(FBTF*LFBTF+FBTA*LFBTA)/FBT 

READ  (VFO    LCG?)LVF() 

PRNT 

PRNT 

PRNT 

PRNT 

PRNT    VARIABLE    LOAD    RECAP: 

PRNT  I  ["EM  WEIGHT  LCC 


[VFO] 

IN  FORWARD  TANKS?) FBTF 
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1 1  / 

!PRNT 

l  Id 

!PRNT    CHEW    &    EFFECTS 

ICE] 

[LCE] 

1  19 

!PRNT   SLBM'S 

[MCJ 

ILMC] 

120 

!PRNT    SANITARY    TANKS 

1ST! 

[LST] 

121 

!?RNT    SUMP    LUBE    OIL 

[Lo] 

(  LLO] 

122 

!PRNT   CLEAN    F.O.     (FIXED) 

[CFO] 

[LCFO] 

123 

IPRNT    ***TOTAL   FIXED*** 

[FL] 

fLFL] 

124 

IPRNT    PROV.    &    STORES 

[PSJ 

[LPS] 

125 

!PRNT    DEVITALIZATION    02 

[02] 

[L02] 

126 

PRNT   TOR PE DOS    FWD    ROOM 

[TF] 

[LTF] 

127 

!PRNT    TORPEDOS    AFT    R(X)M 

CTA] 

[LTA] 

128 

PRNT   TACT.    MISSILES    FWD 

[MF] 

[LMF] 

129 

PRNT   TACT.    MISSILES    AFT 

[MA] 

[LMA] 

130 

PRNT   WRT  TANKS 

[WRT] 

[LWRTJ 

131 

.PRNT    RESERVE    ELECTROLYTE 

[EL] 

[LELJ 

132 

PRNT   TOTAL    FRESH    WATER 

[FW] 

[LFW] 

1  33 

PRNT    RESERVE    LUBE   OIL 

[RLO] 

[LRLO] 

134 

!PRNT    FUEL    BALLAST    FWD 

[FBTF] 

[LFBTF] 

135 

PRNT    FUEL    BALLAST    AFT 

[FBTA] 

[LF3TA] 

136 

.PRNT    VARIABLE    FUEL    OIL 

[VFO] 

[LVFO] 

13/ 

READ(ARE   THESE    VALUES    ACCEPTABLE    <Y/N>?    )&R 

138 

IF(AR.EQ."NO")    GOTO   START 

1  19 

PRNT 

1  40 

AFL=LSUB-LFL 

141 

APS=LSUB-LPS 

142 

A02=LSUB-L()2 

143 

ATM=LSUB-LTM 

1  44 

AWRT=LSUB-LWRT 

145 

AEL=LSUB-LEL 

1  46 

AFW=LSUB-LFW 

147 

ARLO=LSUB-LRLO 

I4d 

AFBT=LSUB-LFBT 

149 

AVF()=LSUB-LVF() 

150 

ATF=LSUB-LTF 

151 

ATA=LSUB-LTA 

152 

AMF=LSUB-LMF 

153! 

AMA=LSUB-LMA 

154 

AFBTF=LSUB-LFBTF 

155! 

AFBTA=LSUB-LFBTA 

156 

MFL=FL*AFL 

157! 

MPS=PS*APS 

158 

M()2=02*A02 

159 

MTM=TM*ATM 

160 

MWRT=WRT*AWRT 

161 

MEL=EL*AEL 

162! 

MFW=FW*AFW 

163 

MRLO=RLO*ARLO 

164! 

MFBT=FBT*AFBT 

165! 

MVF0=VF0*AVF0 

1  60! 

MTF=TF*ATF 

167! 

MTA=TA*ATA 

1  6d! 

MMF=MF*AMF 

169! 

MMA=WA*AMA 

17U! 

«FBTF=FBTF*AFBTF 

1  l\  ! 

«FBTA=FBTA*AFBTA 

172! 

CVL=.9y375 

1  /3! 

CVH=1  .0093  15 

1  74! 

GOTO    NSTM 

1  /b! 

•STAT    CONTINUE 

1  /6! 

HUNT 

•99- 


!  II 
I  76 

I  /O 

lao 

Idl 

182 

183 

184 

185 

186 

187 

I  88 

189 

100 

101 

102 

103 

104 

105 

106 

107 

108 

I  00 

200 

20.1 

202 

203 

204 

205 

206 

20  7 
208 
200 
210 
211 
212 
213 
214 
215 
216 

21  7 
218 
210 
220 
221 
222 
■>2l 
224 
■>2b 
226 
221 
■>2H 
220 
230 
23  1 
232 
233 
234 
235 
236 


PRNT    BALANCE    OF   VARIABLE    LOAD    ACCOUNT*     [BAL] 

PRNT 

RTNSUB 

#NSTM    CONTINUE 

LDSUB=W( |4)-H ( 10) 

#HEAVI  CONTINUE 

CW=4  3/35 

H I =FL+TM+. 5*PS+EL+WRT+. 5*RL0+. 5*FW+CW*FBT 

MH1  =MFL+MTM+.  5*MPS+MEL+MWRT+. 5*MRL0+. 5*MFW+CW*MFBT 

VBHI=LDSUB-MBT-HI 

AMBT=L5UB-LM3T 

MMBT=MBT*AMBT 

MVBHI=-MMBT-MH1 

XV  (I )=MVBHI*CVL 

YV(  I  )=VBHI*CVL 

#HEAV2   CONTINUE 

H2=FL+TM-t-PS+EL+WRT+RL()+FW+CW*FBT 

MH2=MFL+MTM+MPS+MEL+'>1WRT+MRL0+MFW+CW*MFBT 

VBH2=LDSUB-MBT-H2 

MVBH2=-MMBT-MH2 

XV(2)=MVBH2*CVL 

YV(2)=VBH2*CVL 

#LITEI 

LI=FL+.75*PS+WRT+.75*RLO+. /5*FW+FBT+VF0 

MLI=MFL+.75*MPS+MWRT+. /5*MRL0+. /5*MF*+MFBT+MVF0 

VBLI=LDSUB-MBT-LI 

MVBL1=-MMBT-ML1 

XV(3)=MVBLl*CVH 

YV(3)=VBL1*CVH 

#LITE2 

L2=FL+WRT+RL0+VF(HCW*FBT 

ML2=MFL+MWRT+MRL0+MVF()+CW*MFBT 

VBL2=LD8UB-MBT-L2 

MVBL2=-MMBT-ML2 

XV(4)=MVBL2*CVH 

YV(4)=VBL2*CVH 

#HF1 

HF1=FL+MF+TF+.  75*PS*EL+WHT+.  75*RLO+.  75*FW  +  VF()+FBTA+CW*FBTF 

MHFI=MFL+MTF+MMF+.  /5*MPS+MEL+MWRT+.  /5*MRL()+.  /5*MFW+MVFo+MFBTA+CW*WFBTF 

VBHF1 =LDSUB-MBT-HF1 

MVBHFI=-MMBT-MHFI 

XV(5)=MVBHFI 

YV(5)=VBHFI 

#HF2 

HF2=FL+MF+TF+.5*Pd+EL+WRT+.25*RLO+.5*FW+VFO+FBTA+CW*FBTF 

MHF2=MFL+MTF+MMF+.75*MPS+MEL+MWRT+.25*MRLO+.5*MFW+MVFO+MFBTA+CW*MFBTF 

VBHF2=LDSUB-MBT-HF2 

MVBHF2=-MMBT-MHF2 

XV(6)=MVBHF2 

YV(6)=VBHF2 

#HA 

HA=FL+MA+TA+.5*PS+WRT+.  75*RLO+.5*FW+VFO+FBTF+CW*FBTA 

MHA=MFL+MMA+MTA+.5*MPS+MWRT+.75*«RL()+.5*MFW+MVFO+MFBTF+C^*MFBTA 

VBHA=LI)SUB-MBT-HA 

MVBHA=-MMBT-MHA 

XV( 7)=MVBHA 

YV (  /)=VBHA 

#N 

N=H2 

MN=WH2 
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23/ 

238 

239 

240 

24  1 

242 

24  3 

244 

245 

246 

247 

24d 

24Q 

250 

251 

252 

253 

254 

255 

250 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

266 

269 

2/0 

271 

2/2 

273 

274 

2/5 

276 

277 

2  7ti 

2/9 

280 

281 

282 

263 

204 

285 

286 

267 

2  80 

289 

290 

29| 

2Q2 

29  3 

29  4 

295 

2^6 


VBN=LDSUB-MBT-N 

MVBN=-MMBT-MN 

XV(8)=MVBN 

YV(6)=VBN 

#C()NDM   CONTINUE 

CM=H2-CW*FBT+FBT 

MCM=MH2-CW*MFBT+MFBT 

VBCM=LDSUB-MBT-CM 

MVBCM=-MMBT-MCM 

XV(9)=MVBCM 

YV(9)=VBCM 

<#*************  SET  UP  LOCATION  OF  TRIM  &  AUX  TANKS  *************** 

PRNT 

PRNT  YOUR  ACTUAL  DESIGN  MAY  HAVE  AS  MANY  TRIM  AND  AUX  TANKS  AS  YOU  DEEM 

PRNT  APPROPRIATE  FOR  YOUR  PARTICULAR  ARRANGEMENT  SCHEME.  TO  SIMPLIFY  THE 

PRNT  CALCULATION  PROCESS  FOR  THE  EQUILIBRIUM  POLYGON.  YOU  ARE  ASKED  TO 

PRNT  ENTER  AGGREGATE  CAPACITY  AND  LCG  FOR  THE  FORWARD  TRIM  GROUP.  THE 

PRNT  AFTER  TRIM  GROUP.  AND  THE  AUXILIARY  GROUP  TANKS. 

PRNT 

GOTO  SIZE 

#L(X)P 

EXECV 

READ(LCG  OF  FORWARD  TRIM  TANK  GR0UP7JLFT 

READ(LCG  OF  AFTER  TRIM  TANK  GR()UP?)LAT 

READ(LCG  OF  AUX  TANK  GH()UP7)LAUX 

AFT=LSUB-LFT 

AAT=LSUB-LAT 

AAUX=LSUB-LAUX 

RTNSUB 

#SIZE  CONTINUE 

J=l 

REPEAT 

IF( XV(J) .GT.5000)GOTO  LARGE 

IF(XV(J).LT.-5000)G0T0  LARGE 

IF(YV(J).GT.500)G0T0  LARGE 

J=J-H 

UNTILU.GT.9) 

GOTO  SMALL 

#LARGE 

DEL   PAR    MIT.HULLOUT 

EXIT  PART  FILE  MIT.HULLOUT 

XS=IOOO 

YS=100 

ACT  PAR  MIT.POLYOUT 

ACT  DRA  POLYGON  FORM  MIT. LARGEPoLY  DRAW  PoLY 

SEL  CPL  LEFT 

DEF  VIE  GRAPH:  X  I  1 Y6. XOYO, X22Y 17  <CR> 

GOTO  GRAPH 

#3MALL   CONTINUE 

DEL  PAR  MIT.HULLOUT 

EXIT    PART    FILE    MIT.HULLOUT 

XS=500 

YS=50 

del  par  mit.polyout 

act  par  mit.polyout 

act  dra  polygon  form  mit.smallpoly  draw  poly? 

sel  cpl  left 

def  vie  graph*  x 1 1 y6.x0y0. x ?2y i  /  <cr> 

#graph  continue 

J=l 
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297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

31b 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

130 

33  1 

332 

333 

334 

135 

336 

33  7 

33a 

339 
340 
341 
342 
343 
3  44 
345 
346 
347 
343 
349 
350 
351 
3i>2 
353 
354 
3  55 
356 


REPEAT 

XP(J)=XV(J)/XS 

YP(J)=YV(J)/YS 

INSERT  POINT'  XIXP(J) JYrYP(J)]  <CR> 

XT(J)=XP(J)+.2 

YTU)=YP(J)-.  I 

J=J  +  I 

UNTIL(J.GT.9) 

INS  TEX  HI  HG  .2«  X[ XT < I ) I YI YT ( I ) )     <CR> 

INS  TEX  H2  HG  .2«  XI XT(2) I Yl YT(2 ) J  <CR> 

INS  TEX  LI  HG  ,2»  XI  XT (3) ] YI YT (3 ) J  <CR> 

INS  TEX  L2  HG  .2:  XIXT( 4) ] YI YT< 4 ) ]  <CR> 

INS  TEX  HFI  HG  .2«  X[ XT(5) JYf YT( 5) J  <CR> 

INS  TEX  HF2  HG  .21  X[XT(6 ) JY[ YT( 6 ) ]  <CR> 

INS  TEX  HA  HG  .2«  XI XT( 7) ] Yl YT( 7 > ]  <CR> 

INS  TEX  N  HG  .2:  XI  XT <d )-.5 ] YI YT (a ) ]  <CR> 

ins  tex  m  hg  .2«  xi xt(9 ) ] yi yt(9 ) ]  <cr> 
gosub  Loop 

GOTO  PLTNK 

#CAP  CONTINUE 

EXECV 

REALMCAPACITY  <TONS>  OF  FORWARD  TRIM  GROUP?)CFT 

REAO(CAPACITY  <TONS>  OF  AUXILIARY  GR()UP?)CAUX 

READ(CAPACITY  <T()NS>  OF  AFTER  TRIM  GROUP?)CAT 

RTNSUB 

#PLTNK    CONTINUE 

GOSUB   CAP 

M  (  I  ) =0 

M(2)=AFT*CFT 

M(3)=M(2)+AAUX*CAUX 

M(4)=M(3)+AAT*CAT 

M(5)=M(4)-AFT*CFT 

M(6)=M(5)-AAUX*CAUX 

M ( 1 ) =0 

Y( I )=0 

Y(2)=CFT 

Y(3)=CAUX+CFT 

Y(4)=Y(3)+CAT 

Y(5)=Y(4)-CFT 

Y(6)=Y(5)-CAUX 

Y  (  / )  =0 

J=l 

REPEAT 

M(J)=M( J)/XS 

Y(J)=Y( J)/YS 

J=J+I 

UNTIL(  J.GT.7) 

INS    LIN   TAG=I«    XfMd  )  JYIYd  )  ].XfM(2)  )YIY(2)  1  ,    <CR> 

INS    LIN   TAG=2:    XI M( 2) ] YI Y( 2 ) 1 . XI M( 3 ) J YI Y( 3 ) ] .    <CR> 

INS    LIN   TAG=3«     XTM (3 ) 1 YI Y (3) J , XI M (4 ) ] YI Y (4 ) ] ,     <CR> 

INS    LIN   TAG=4«    X  t  M  <  4  )  ]  Yl  Y  (  4  )  ]  .X I  M<  5  >  ]  Yl  Y(  5  )  J  ,    <CR> 

INS    LIN   TAG=5*    X f M (5) I Yt Y (5 ) 1 . X[M (6 ) ] Yt Y (6 ) ] .     <CR> 

INS    LIN    TAG=6:    X[M( ft ) J YI Y( 6) J . XI M( 7)  J Y[ Y( 7  )  J  .    <CR> 

EXECV 

READ*  IS    THIS    POLYGON    SATISFACTORY    <Y/N>       ?       )AR 

IFt&W.EQ.  "Y">    GOTf)    FINAL 

DEL    ENT*    TAG    I.TAJ    2. TAG    3. TAG    4, TAG    5. TAG   6    <CR> 

EXECV 

REAIXCHANGE  CAPACITIES  ONLY  <l>.  OR  LCG'S  &    CAPACITIES  <2>  ?  )CHC 

IF(CHC.H0.2)  GOTO  BOTH 
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357 
.158 
359 
300 
361 
362 
363 
364 
305 
366 
367 
368 
369 
3/0 
371 
3/2 
373 
374 
375 
376 
37/ 
378 
3/0 
380 
381 
382 
3d3 
384 
3d5 
386 
387 

3  88 
369 
390 
39  1 
392 
393 
394 
395 
396 
39/ 
398 
3QQ 
400 
401 
402 
403 
404 
405 
406 
407 
406 
409 

4  10 
4  I  I 
412 
413 
414 
415 
416 


IF(CHC.EQ. I )    GOTO    PLTNK 

#BOTH    CONTINUE 

GOSUB    LOOP 

GOTO    PLTNK 

#FINAL   CONTINUE 

OPENW    3.J'MIT.PASS.P0LYDAT" 

&X=CE 

GOSUB    O 

&X=LCE 

GOSUB    O 

&X=MC 

GOSUB    0 

&X=LMC 

GOSUB    O 

&X=ST 

GOSUB    0 

&X=LST 

GOSUB    0 

&X=LO 

GOSUB    0 

AX=LL<> 

GOSUB    0 

&X=CF() 

GOSUB    0 

&X=LCFO 

GOSUB    0 

&X=FL 

GOSUB    0 

&X=LFL 

GOSUB    O 

&X=PS 

GOSUB    0 

&X=LPS 

GOSUB    O 

&X=02 

GOSUB    O 

AX=L02 

GOSUB    0 

&X=TF 

GOSUB    0 

&X=LTF 

GOSUB    0 

&X=TA 

GOSUB    0 

&X=LTA 

GOSUB    O 

&X=MF 

GOSUB    O 

&X=LMF 

GOSUB    0 

&X=MA 

GOSUB    0 

&X=LMA 

GOSUB    O 

*X=WHT 

GOSUB    0 

&X=LWWT 

GOSUB    0 

&X  =  EL 

30SUB    0 
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4|7!&X=LEL 
41b!G0SUB  0 
4|P!&X=FW 
42U1G0SUB  O 
421 !&X=LFW 
4  22SGOSUB  0 
423!&X=KLO 
4241GOSUB  0 
425!&X=LHL() 
426!G0SUB  0 
427!&X=FBTF 
428SGOSUB  0 
4  29!&X=LF3TF 
4301GOSUB  O 
431  !&X=FBTA 
432SGOSUB  O 
433!&X=LFBTA 
4341G0SUB  0 
435!&X=VFO 
436IGOSUB  0 
437!&X=LVFO 
43d!G0SUB  <> 
439!&X=CFT 
440!G0SUB  0 

44  1 !&X=LFT 

442IG0SUB  0 

443!&X=CAUX 

444!G()SUB  0 

445!&X=LAUX 

446SG0SUB  t) 

447!&X=CAT 

44d!G0SUB  0 

44Q!&X=LAT 

4501G0SUB  O 

451 !PLOT  DOT  SCALE  .5 

4521EXIT  PART  F 

453!G()TO  END 

454!  #<) 

455!WRITEF  3.&X 

4561RTN5UB 

457!#END 

4581RUN  NEW  MIT.POLYFILfc 

459! END 
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I 

2 

3 

4 

5 

6 

7 

8 

Q 

10 

I  I 

12 

13 

14 

15 

16 

1  7 

Id 

10 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

30 

40 

41 

42 

43 

44 

45 

46 

4  / 

48 

40 

50 

■j\ 

52 

-j  3 

54 

5b 

56 


OPENR  I ,»MIT. PASS.POLYDAT" 

OPENW  2. "MIT. OUT. POLY" 

PRNT  OUTPUT  FOR  THE  POLYGON  MODULE  WILL  BE  IN  'MIT. OUT. POLY', 

AT="P0LYG0N  SUMMARY  OUTPUT  FILE" 

WRITEF  2. AT 

AP=»       '• 

WRITEF  2.AP 

WRITEF  2,&P 

AT="     ITEM  LCG  WEIGHT" 

WRITEF  2, AT 

WRITEF  2,AP 

READF  1.AX 

READF  1  ,AY 

AT="CREW  ANO  EFFECTS  '<+&Y+J'  "+&X 

WRITEF  2. AT 

GOSUB  RD 

&T="SLBM'S  OR  COMP.  WATER      "+AY+"  "+&X 

WRITEF  2, AT 

GOSUB  RD 

&T="SANITARY  TANKS  "+AY+"  "+AX 

WRITEF  2, AT 

GOSUB  RD 

AT="LUBE  OIL  IN  SUMPS  "+&Y+"  J'+&X 

WRITEF  2. AT 

GOSUB  RD 

AT="FIXED  CLEAN  FUEL  OIL       »+aY+"  "  +  AX 

WRITEF  2, AT 

GOSUB  RD 

AT="**** TOTAL  Fi XED  ITEMS****  "+AY+"  "+&X 

WRITEF  2.&T 

GOSUB  RD 

AT  =  "PROVISI')NS  AND  STORES      «+&Y  +  J'  "+&X 

WRITEF  2, AT 

GOSUB  RD 

AT="REVITALIZATION  02  »+&Y+"  "+&X 

WRITEF  2, AT 

GOSUB  RD 

&T="TORPED()S  IN  FWD  ROOM       "+AY  +  "  "+&X 

WRITEF  2, AT 

GOSUB  RD 

AT="TORPEDOS  IN  AFT  ROOM       "+&Y+"  "+AX 

WRITEF  2. AT 

GOSUB  RD 

AT="TACTICAL  MISSILES  FWD      "+AY+"  "  +  AX 

WRITEF  2, AT 

GOSUB  RD 

AT="TACTICAL  MISSILES  AFT      "+AY+"  "+AX 

WRITEF  2. AT 

GOSUB  HI) 

&T="WHT  TANKS  »+&Y+i'  "+&X 

WRITEF  2, AT 

GOSUB  RD 

AT="RESERVE  ELECTROLYTE        "+KY+"  "t-AX 

WW1TEF  2.*T 

GOSUB  RD 

AT="T0TAL  FRESH  WATER  "+XY+"  "  +  KX 
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5/SWRITEF    2. AT 

58 SGOSUB    RD 

59! AT=" RESERVE    LUBE    OIL 

"+AY+" 

60SWRITEF   2. AT 

6J  SGOSUB    RD 

62!AT="FUEL   BALLAST  FWD 

"+AY+" 

63SWRITEF    2,&T 

64 SGOSUB    RD 

65!&T="FUEL    BALLAST   AFT 

"+AY+" 

66!HRITEF    2.&T 

67 SGOSUB    RD 

o8!AT="VARlABLE   FUEL    OIL 

"+AY+" 

6QSWRITEF   2, AT 

70 SGOSUB    RD 

/I SHRITEF    2.AP 

72!rtRITEF    2.&P 

/3!&T="VARIABLE    BALLAST   TANK 

DATA»" 

74SWRITEF    2, AT 

75!&T="         TANK    GROUP 

LCG 

76SWRITEF    2  ,AT 

77!&T="F()RWARD   TRIM 

"+AY+" 

78SWRITEF   2. AT 

79 SGOSUB    RD 

dO!AT="AUXILlARY 

"+AY+" 

8) SWRITEF    2, AT 

B2SG0SUB    RD 

83!&T="AFTER   TRIM 

"+AY+-" 

84SHRITEF    2, AT 

85! GOTO    END 

86!#RU    CONTINUE 

87SREADF    1 ,AX 

88SREADF     1 .AY 

89SRTNSUB 

90!#END 

91  SEND 

"  +  AX 
"+AX 
"+AX 
"+AX 


CAPACITY    (TONS)" 
"+AX 

"+AX 

"+AX 
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